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^ . Abstract: We analyze the renormalization of systems whose effective degrees of free- 

^ ' dom are described in terms of fluctuations which are "environment" dependent. Relevant 

environmental parameters considered are: temperature, system size, boundary conditions, 
and external flelds. The points in the space of "coupling constants" at which such systems 
exhibit scale invariance coincide only with the flxed points of a global renormalization 
I group which is necessarily environment dependent. Using such a renormalization group we 

give formal expressions to two loops for effective critical exponents for a generic crossover 
induced by a relevant mass scale g. These effective exponents are seen to obey scaling laws 
^ ■ across the entire crossover, including hyperscaling, but in terms of an effective dimension- 

ality, d^ff — 4 — 7;^, which represents the effects of the leading irrelevant operator. We 
analyze the crossover of an 0{N) model on ad dimensional layered geometry with periodic, 
antiperiodic and Dirichlet boundary conditions. Explicit results to two loops for effective 
exponents are obtained using a [2,1] Pade resummed coupling, for: the "Gaussian model" 
[N = —2), spherical model (A^ = oo), Ising Model {N = 1), polymers [N = 0), XY-model 
{N = 2) and Heisenberg [N = 3) models in four dimensions. We also give two loop Pade 
resummed results for a three dimensional Ising ferromagnet in a transverse magnetic fleld 
and corresponding one loop results for the two dimensional model. One loop results are 
also presented for a three dimensional layered Ising model with Dirichlet and antiperiodic 
boundary conditions. Asymptotically the effective exponents are in excellent agreement 
with known results. 



§ 1. INTRODUCTION 



Field theory is our most powerful tool for describing systems with a large number of 
degrees of freedom, and in particular in situations where the latter can act collectively. In 
one guise or another such systems comprise a large proportion of current physics. Even 
though one starts at a certain energy scale with a "fundamental" set of degrees of freedom, 
which may be very elementary looking, e.g. QED, one knows that as one progresses to 
different (usually lower) energy scales, the original elementary degrees of freedom, though 
in principle offering a complete description of the low energy physics, become, in practice, 
unusable. 

For example, one would be rather perverse to try and describe a phase transition 
in a ferromagnet in terms of QED. As is well known the transition is much better de- 
scribed in terms of a Landau- Ginzburg- Wilson Hamiltonian (LGW Hamiltonian) based 
on effective degrees of freedom associated with a magnetization density (f){x). This is an 
implementation of the philosophy of effective field theory, whereby the theory provides a 
good description only at energy scales where the effective degrees of freedom chosen of- 
fer a reasonably faithful representation of the physics. For instance, if we heated up the 
ferromagnet until it became a plasma of nuclei and electrons, then a description in terms 
of a LGW effective Hamiltonian would be rather inappropriate. The difficulty, of course, 
is that the effective degrees of freedom of the system are scale dependent, i.e. there is a 
"crossover" in them, and consequently in the physical behaviour of the system. Thus, if 
one requires a complete description of the physics at all scales one must be able to account 
for this fact. Our aim in this paper will be to develop further a formalism for describing 
such systems using field theoretic renormalization groups (RG) . We will for the most part 
keep the formalism as general as possible, specializing in the latter half of the paper to 
specific crossovers of interest. Our reason for doing this is that the formalism is applicable 
to a very wide class of crossovers. By presenting the general results the reader is then at 
liberty to treat a crossover of particular interest to him/her merely by inserting into the 
appropriate parts of the formalism the particular details specific to that crossover. We will 
use a language that is hopefully accessible to people with some field theoretic background 
in both particle physics and statistical physics as the applications of the formalism are 
equally applicable to both 

To illustrate the ubiquity of the concept of crossover behaviour we will briefiy men- 
tion (without detailed explanation) some physically pertinent examples. In QCD, at high 
energies, a description in terms of quark-gluon degrees of freedom is appropriate, whilst 
at low energies baryon-meson degrees of freedom are more suitable. In four dimensional 
finite temperature field theory, for T ^ m{T), where m{T) is a typical finite temperature 
mass scale, three dimensional degrees of freedom are appropriate, whereas for T <^ m(T) 
they are four dimensional. In the early universe the effective degrees of freedom are time 
dependent due to the cosmological expansion. In fact, for quantum field theory in curved 
spacetime where one is concerned with "in" and "out" vacua, the difference between these 
vacua, which results in particle production, can be thought of as being representative of 
the difference between "in" and "out" effective degrees of freedom. In critical phenomena 
crossover can be induced by many different effects. One exciting much interest currently 
is the effect of randomness, e.g. a random magnetic field, random impurities etc. In such 
cases the effective degrees of freedom depend on the concentration of the impurities or the 
strength of the magnetic field. Another eliciting great interest concerns the effects of finite 
size. This is of relevance in lattice simulations and also in real experimental systems. 

The physics of crossovers then, concerns systems which exhibit qualitatively different 
degrees of freedom at different scales. Our goal is to qualitatively and quantitatively 
describe such systems. For an exact model of course there is no problem, however, even 
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though one has an exact answer that does not necessarily imply that one has a good 
intuitive understanding of the physics. Often with an exact model one can not see the wood 
from the trees. Generically one must resort to approximation techniques. The dominant 
one, of course, is perturbation theory, where one must answer the crucial question: what 
parameter should one perturb in? Often the difference between getting a nonsensical 
versus a reliable answer comes down to being able to pick and work with a good expansion 
parameter. 

In crossover problems one starts with some interacting Hamiltonian which will come 
with a coupling constant, A, which one would be tempted to use as an expansion parameter 
when computing the correlation functions associated with this Hamiltonian. However, the 
starting "bare" (microscopic) Hamiltonian will only offer an accurate perturbative descrip- 
tion when one looks at scales k ~ A, A being the UV cutoff (inverse lattice spacing) for 
the theory. At scales k ^ A, as is well known, the bare parameters offer a perturbatively 
useless description. One must therefore perform a renormalization to a new set of param- 
eters. This is done by defining correlation functions at a fiducial scale k', and using them 
as the parameters with respect to which one describes physics at a scale k ^ k' . As long 
as the two scales k and k' are not too different, the description of physics at the scale k in 
terms of parameters defined at the scale k/ will be perturbatively reliable. The crossover 
from bare degrees of freedom to renormalized ones, associated with taking the continuum 
limit or the removal of a regularizing UV cutoff, is perhaps the most elementary example 
of a crossover. In contrast to a scale invariant system where the bare degrees of freedom 
and the renormalized ones are relatively simply related they are quite different when the 
system is not scale invariant. In the presence of a cutoff the bare and renormalized degrees 
of freedom are always different as the system cannot be scale invariant. 

So why are the bare and renormalized effective degrees of freedom so different? Due 
to the effects of fluctuations; the fluctuations "dress" the bare parameters. This dressing 
in principle can become infinite if one takes the A — (X) limit. For us the crucial difference 
is between a large dressing and a small dressing not between an infinite one and a finite 
one. Moreover, as far as the RG is concerned, we feel that there is no particular virtue 
in emphasizing UV dressings as opposed to IR ones. In particle physics it is the former 
which have received most attention, however, the machinery of the RG must be invoked 
any time one has large dressings of one's parameters, irrespective of from which particular 
end of the spectrum the dominant fluctuations arise. Large dressings imply that one is 
not tracking the effective degrees of freedom in the problem very well. Renormalization 
is a methodology by which one can use the freedom to reparametrize one's theory to find 
parameters which track the effective degrees of freedom more accurately. In the case of 
bare versus renormalized parameters, without interactions the bare parameters cannot get 
dressed and remain the same. Indeed, it is the presence of interactions that induces the 
crossover between them. The renormalization carried out was dependent on the interaction 
strength A. That is why the renormalization worked in the first place, because it depended 
on the parameter that was inducing the crossover. We could have been perverse and tried 
to renormalize the sytem in a A independent way. The Gaussian theory would have a large 
dressing associated with the specific heat, or vacuum energy, for instance. This would have 
to be renormalized. After the rcnormalizations associated with the Gaussian theory were 
performed one would still find that connected correlation functions had large dressings. 
This is a refiection of the fact that the renormalization process, using only counterterms 
appropriate for the Gaussian theory, was not sufficient to make sense of the interacting 
theory. 

Generically, in the renormalization process one takes fiuctuations between scales A and 
contributing to a given correlation function, and absorbs them into a redefinition of 
a parameter, which is usually a vertex function at the scale k! . To investigate physics 
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at a scale k, one considers the fluctuations between k and k' , explicitly in perturbation 
theory, in terms of the coupling A defined at the scale k' instead of the bare coupling. The 
important question of course is what fluctuations are being absorbed into redefinitions of 
the parameters? If we renormalized using only counterterms associated with the Gaussian 
theory then only non-interacting fluctuations are being absorbed into redefinitions. For 
the connected correlation functions this is equivalent to not having absorbed in anything. 
Obviously this won't help much. In the standard A dependent case all the fluctuations 
associated with interacting effective degrees of freedom between the scales A and k' have 
been absorbed. What we do know for a fact is that the effective degrees of freedom 
in the system are interacting, therefore one would expect to get more sense out of a 
renormalization procedure that explicitly takes account of this fact. The reader might 
wonder why we have gone to such lengths analyzing such a standard procedure, and why 
on earth one would ever try to renormalize an interacting theory using renormalization 
that was explicitly A independent. We hope by the end of the paper the reason will be 
clear. For the moment we will simply say that a successful renormalization is one that can 
take proper account of what the true effective degrees of freedom of a system are. In the 
above they are A dependent, therefore a good renormalization should be A dependent. 

So, an important aspect of renormalization is that it is a methodology whereby through 
a suitable redefinition of parameters one can make perturbative sense out of something that 
originally was perturbative nonsense. How sensible the expansion actually becomes, how- 
ever, is crucially dependent on what physics the renormalization process can capture. In 
the above the renormalization had to capture the fact that the effective degrees of free- 
dom were interacting. If one considered an interacting field theory in a three dimensional 
box of size L, one could renormalize the theory in an L independent fashion. When one 
considered physics on scales k ~ one would find that the theory was perturbatively 
ill defined, whereas an appropriate L dependent renormalization made perturbative sense. 
The reason for this, of course, is that the effective degrees of freedom in the system are 
explicitly L dependent. An L independent renormalization ignores this important physical 
fact. The only fluctuations being absorbed into the renormalized parameters in this case 
are L independent, no matter what renormalization scale one chooses. L here is the pa- 
rameter which induces the crossover and therefore a good renormalization scheme should 
be L dependent. 

More generally we will speak of a crossover being induced by a particular "environmen- 
tal" variable. The idea to understand here is that the fluctuations in a system "feel" out 
the environment they are in. For instance in the above case of physics in a box, the size 
of the box L is a relevant environmental variable, as clearly effective degrees of freedom 
associated with length scales <^ L and ~ L are completely different. The effective degrees 
of freedom arc thus sensitive to the environment. There are many, many different environ- 
mental parameters describing all sorts of crossover systems. In this paper we will consider 
the effects of a generic environmental variable g, which could represent variously: finite 
size, temperature, spin anisotropy, long range interactions, cosmological constant, Hubble 
parameter, magnetization, electric/magnetic field and so on. What we will be showing is 
how to describe qualitatively and quantitatively the response of a system to changes in 
the environment. We will show that in order to obtain a perturbative description of such 
systems one must implement an RG which is explicitly dependent on the environment. 
We will call such RGs "environmentally friendly". If an environment "unfriendly", that 
is to say, environment independent, RG is implemented, we will see that this generically 
leads to breakdown of perturbation theory, wherein the theory becomes strongly coupled 
in terms of the effective degrees of freedom characteristic of the environmentally unfriendly 
RG, and that large perturbative dressings appear which can actually become divergent. 
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Crossover behaviour, as mentioned, is quite ubiquitous. It is not our intention here 
to give a review of what has been done, consequently we will be rather selective in our 
comments, and, more apologetically, in our references. In critical phenomena it has recieved 
much more attention than in particle physics (for an early review see [1]), though many 
problems in particle physics are essentially crossover problems. In principle, basically any 
laboratory system will exhibit crossover behaviour in some regime. Some of the more 
experimentally accessible ones are: uniaxial dipolar ferromagnets [2], systems exhibiting 

a bicritical point [3], dimensional crossover in liquid He^ [4], quantum ferromagnets ([5] 
and references therein), bulk/surface crossovers [6]. Dimensional crossover, one of the chief 
concerns of this paper, has been studied mainly in the context of finite size scaling [7]. 
Much work has been done on the latter, and finite size effects in general, in the context of 
lattice simulations; quantum Monte Carlo methods [8] are also closely related. 

From a more theoretical standpoint, and more particularly from the point of view of 
RG theory, much work has also been done. For a typical crossover induced by a generic 
anisotropy parameter gr, scaling formulations have been almost invariably based on RG's 
that are g independent [9]. This automatically leads to critical and crossover exponents 
that are defined with respect to the isotropic g = Q fixed point. Such RG's would be 
incapable, in and of themselves, of bridging the crossover to the anisotropic g = oo fixed 
point. What we mean by this is that the fixed points of such RG's will not yield all 
the possible points of scale invariance of the system, in particular the anisotropic scale 
invariant system will be inaccessible. The fundamental length scale utilized is ^q, the 
correlation length of the isotropic system. It is pertinent to note that this is not the physical 
correlation length in the crossover system. The desire to make accessible another fixed 
point besides the isotropic one has often entailed the matching of asymptotic expansions 
around the anisotropic and isotropic fixed points, or the use of high temperature expansions 
in conjunction with an ansatz for the scaling function [10] [11]- Exact models have also 
played a role, for instance, the two dimensional Ising model [12] and the spherical model 
[13] in the context of finite size scaling. 

An RG approach, due to Riedel and Wegner [14], and used by others, utilizes "model 
recursion relations" . Here it is assumed known what the two fixed points of the crossover 
are, subsequently RG equations for the scaling fields are postulated that interpolate be- 
tween the two fixed points. As there is no reference to the underlying microscopic Hamilto- 
nian such a ploy can only offer qualitative information. In addition there are many different 
crossover systems that exhibit crossover between exactly the same two fixed points. Uni- 
versality implies that two systems that lie in the same universality class can via a suitable 
rescaling of variables be shown to have exactly the same IR properties, e.g. a simple fluid 
and an Ising ferromagnet. For crossovers the asymptotic fixed points do not determine 
"crossover universality classes" . The scaling functions associated with different crossovers 
between the same two fixed points cannot generally be transformed into one another by a 
change of scaling variable. However, there is a concept of two systems being in the same 
crossover universality class — a (d — l)-dimensional Ising model in a transverse magnetic 
field and a d-dimensional layered Ising model with periodic boundary conditions being a 
case in point. Wilsonian RGs of the "momentum shell integration" type have been fre- 
quently employed with varying success. These methods generally have two main defects: 
ease of extension beyond first order, and using too rough an approximation in the shell 
integration. Given that many physically different systems can crossover between the same 
asymptotic fixed points it is important to be able to pick out the details of the crossover 
curves in order to distinguish between them. By using an approximation on the momentum 
shell integration it is easy to blur such distinctions. 

From a field theory point of view much less work has been done. This is mainly because 
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of the preoccupation with renormahzation in the context of UV divergences, though there 
are one or two notable exceptions. Amit and Goldschmidt [15] introduced the concept of 
generalized minimal subtraction (GMS) in the context of crossover at a bi critical point. 
Their results for 7e/j, however, differ significantly from those found by Nelson and Domany 
[16] and independently by Seglar and Fisher [17] using momentum shell integration, in 
particular in that the latter find a characteristic "dip" in the curves. In this case the 
results of our methodology applied to a bicritical crossover [18] agree with the latter. We 
believe this to be due to a deficiency of GMS which fails to capture the true behaviour of 
the leading irrelevant operator which in our crossover formalism plays a crucial role. GMS 
was also applied to uniaxial dipolar ferromagnets in [19], once again the results of our 
analysis [20] [18] are somewhat different. Schmeltzer [21] calculated 7e// to one loop for 
three dimensional quantum ferroelectrics. Lawrie [22] considered dimensional crossover for 
cZ-dimensional quantal and d+ 1 dimensional finite-sized Ising models for 3 < (i < 4. Unlike 
our methods the s expansion methodology he used could not capture the crossover between 
two non-trivial fixed points as such an expansion is around the upper critical dimension 
which changes across the crossover. Nemirovsky and Freed ([23] and references therein) 
used minimal subtraction techniques but failed to be able to access the full crossover. Field 
theoretic results for dimensional crossover in a fully finite geometry or a cylinder have been 
obtained [24] but the techniques used have not been extended to the case of a system with 
more than one fixed point. In a particle physics context some interesting work has been 
done in the context of the finite temperature RG [25] (see [26] for some recent applications 
of our work to finite temperature field theory). 

In this paper we will present a field theoretic formalism which is applicable to a very 
wide class of crossovers associated with a field theory in a particular environment. We have 
previously presented some one loop results in the case of dimensional crossover, where the 
environmental variable was L, the finite size of the system, both in the case of above [27] 
[28] and below [29] the critical point. In the latter the background magnetization was 
also a relevant environmental variable inducing a crossover. We have also considered, at 
one loop, uniaxial dipolar ferromagnets and bicritical systems in [20] and [18], where the 
environmental variables were ctg ^-nd m respectively. Here is the dipole-dipole coupling 
strength and m is the mass of the 0{N — M) components for the case 0{N) — > 0{M). 
The main thrust of the present paper is to generalize the techniques used to a generic 
environment induced crossover, and in particular, to extend the methods beyond one loop 
as there are many subtleties involved in doing so, not least of which is the application 
of resummation techniques. We present two loop expressions for effective exponents and 
some scaling functions, which are set up in just such a fashion that the reader interested 
in a specific crossover, may take them, insert the relevant environment dependence in 
the Feynman diagrams, crank the handle (analytically if possible, numerically if not), 
and spit out the answers. To explicitly demonstrate the efficacy of the methodology we 
consider dimensional crossover for a d dimensional 0{N) model on a layered geometry 
with periodic, Dirichlet or antipcriodic boundary conditions. We give explicit, surprisingly 
simple, expressions for two loop effective exponents for a four dimensional layered system 
with periodic boundary conditions, and corresponding (not so simple) expressions for a 
three dimensional layered system. We also give results for Dirichlet and periodic boundary 
conditions. Two loop results for the effective exponents of a three dimensional Ising model 
in a transverse magnetic field are presented, as are one loop results for the corresponding 
two dimensional model. Graphical displays of our results can be found throughout the 
paper. 

The format of the paper will be as follows: in section 2 we consider some formal as- 
pects of renormahzation, mainly to set notation. In section 3 we consider the relationship 
between the RG and the environment. Specifically, we compare and contrast the field 
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theoretic and "Wilsonian" RGs asking and answering the question of what does one re- 
quire of a "good" RG, and in particular addressing which RGs will access all the points 
of scale invariance of a physical system. We illustrate these ideas by comparing massless 
versus massive renormalization showing that the former yields an RG which has fewer fixed 
points than the latter. We discuss precisely why a good RG should be environmentally 
friendly. In sections 4 we derive formal two loop expressions for the Wilson functions of 
our generic g dependent crossover. In section 5 we discuss further the idea of an environ- 
mentally friendly versus environmentally unfriendly RG, showing that the former offers 
a more global description of the crossover sytem, the latter being incapable of accessing 
some of the points of scale invariance of the system. We also discuss an RG where the 
environment itself "runs", which is useful for addressing questions such as: what is the 
shift in critical temperature due to changes in the environment? In section 6 we consider 
the consequences of an environmentally friendly RG, defining effective critical exponents 
and showing that they obey scaling laws, including eflFective hyperscaling with respect to 
an effective dimensionality which is a measure of the importance of the leading irrelevant 
operator across the crossover. We exhibit formal scaling forms for the crossover equation 
of state and crossover coexistence curve in terms of effective exponents. We exhibit non- 
linear scaling fields which interpolate across the crossover becoming at the endpoints of 
the crossover the associated linear scaling fields. We discuss the concepts of crossover uni- 
versality and crossover universality class. In section 7 we address the important question 
of "what should we perturb in?" We argue that an appropriate perturbation parameter is 
the resummed solution of the fi function of the environmentally friendly RG. We empha- 
size that in the context of the RG, perturbation theory is carried out at the level of the 
Wilson functions, all of which as perturbative expansions should be resummed. We resum 
these series using a [2,1] Fade approximant which is appropriate at the two loop level. 
When the Wilson functions are integrated (exponentiated) the resulting functions should 
not then be themselves expanded. We show that an environmentally friendly RG is most 
potent when all environment dependent fiuctuations are absorbed into reparametrizations 
of the parameters. In section 8 we discuss the floating fixed point which yields a good 
approximation to crossovers via the solution of the crossover (3 function as an algebraic 
equation instead of as a differential equation. In section 9 we present two loop diagramatic 
representations for effective exponents and scaling fields. These can be treated as "black 
box" type expressions where one inserts the particular environmental dependence in the 
Feynman diagrams and reads off the corresponding exponent. In section 10 we consider 
one and two loop results for an 0{N) model on a layered geometry with periodic boundary 
conditions paying special attention to some models of interest. We also discuss the effects 
of boundary conditions on the crossover. In section 11 we consider the quantum/classical 
crossover mad dimensional Ising model in a transverse magnetic field. We present both 
one and two loop results. In section 12 we draw some conclusions and make some specu- 
lations. There are two appendices containing various expressions for diagrams needed in 
the calculations. 



§ 2. FORMAL RENORMALIZATION 

In this section we present a formal approach to the problem of renormalization. We begin 
with the action (LGW Hamiltonian) 



S[^b] = J d'^x 



111 A 

-{dLpsf + -^rn\^p\ + -tB{x)ip% + -^^b " Hb{x)^B + 9bOb 



(2.1) 



which represents the "microscopic" theory. The term qbOb is used to abstractly refer to 
an "anisotropy" in the system, Ob being an operator conjugate to the coupling qb- We 



7 



will think of (/^ = as representing an isotropic system. Physically this anisotropy could 
be of diverse origin. The general methodology espoused in this paper is not confined to 
sytems described by the Hamiltonian (2.1), though the latter is general enough for most 
of our considerations. From time to time, however, we will extend our attention to a 
general Hamiltonian who's parameters live on an abstract "space of coupling constants" 
Ji4. Each of the parameters in (2.1) is associated with a coordinate direction in this space. 
We will for convenience tend to use the terminology associated with viewing (2.1) as the 
Hamiltonian of a magnetic system. 

The partition function Z is defined by the functional integral 



Z = ^-W[JB+HB,tB+ml] ^ J [(^^]e-'^bB]+ / d'^xJB{x)(pB{x) ^2.2) 

where W is the generating functional of connected Greens functions (free energy for pre- 
scribed sources). Jg is to be understood here as some formal source, differentiation with 
respect to which generates the connected Greens functions (correlation functions), and 
which will after generation of the correlation functions be put to zero. will represent 
a real magnetic field in the critical phenomena context. The efi"ective action (free energy 
for prescribed magnetization), r[0^, + m^], is the generating functional of one particle 
irreducible A'"-point vertex functions and is given by the Legendre transform 

T[(I>B, tB + ml] = W[Jb + Hb, tB] + J d^x{JB{x) + Hb{x))<Pb{x) (2.3) 

where c^b — < '^B >j,h— 4>B'^ < VB >Ji 'Pb being the magnetization when J = 0. The 
vertex functions are then defined in terms of functional derivatives oiT[4>B^tB + ml] with 
respect to 0^ and tB- Specifically 



T[cf>B,tB + ml]= J2 ^^f\^B,tB + ml)icf>B-^Bf 



N=0 

oo ^ 

M=0 

where implicit multiple integrations over position or momenta are understood. 

Note that (pB = above the critical temperature when Hb = 0, and hence the quan- 
tities r'^-^) are independent of (pB in this domain. For Hb ^ 0, or below the critical 
temperature, there is a non-zero spontaneous magnetization ^b- Deviations from this 
magnetization arc responses to the externally applied source Jb, they again will be set to 
zero once the vertex functions are obtained. The equation of state is given by 

rW = Hb (2.4) 

when Jb = 0. This is the basic equation of the effective action and serves to determine 

the background magnetization (pB- 

The calculation of correlation functions more often than not requires the introduction of 
an ultraviolet (UV) cutoff (physical or otherwise) . An infrared (IR) cutoff is frequently also 
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necessary. We renormalize by imposing normalization conditions, though for the moment 
we will be somewhat nebulous about the precise normalization point. Two quantities 

required in this procedure are the wave function renormalization constant Z^, where (j) = 

_ 1 

'^(pB, and the composite operator renormalization constant Z^2. The relation between 
the bare and renormalized vertex functions is 



{N,M) ^ ^-f ^-Mj.{N,M) 



B 



(2.5) 



our convention being that the unsubscripted quantity is the renormalized one. The pa- 
rameters of the theory and are multiplicatively renormalized, the relation between 
the bare and renormalized quantities being 



(2.6) 



must be additively renormalized, = + 5rn^. 

The renormalization factors are determined by the following generic normalization 
conditions 



r(2) 



p2 



ar(2) 



p2 



r(2'i) =1 

^ \bP=ai,T=a2 ^ 



(2.7) 
(2.8) 
(2.9) 



r(4) 



SP=ai,T=a2 



A 



(2.10) 



where r = and SP denotes pj^ ■ pj — ^{AS-ij — 1), pi being the momentum entering 
through the jth external leg, and k is an arbitrary renormalization scale. For the moment 
we will not specify the values of the various parameters involved in the normalization 
conditions. The alert reader might wonder why there is no condition corresponding to a 
renormalization of g^- The answer is that in this paper we take it to be a RG invariant 
(i.e. a non-linear scaling field in the sense of Wegner [30]), though this in no way restricts 
the generality of our approach. We define the correlation length via the relation 



2dfddxG('^){x,0) 



(2.11) 



where 



6J{x)dJ{0) 



J=0 



. The conditions (2.7) and (2.8), together with ai = 0, 



imply that ^ defined by (2.11) is t 2 at the normalization point. 

The wavefunction renormalization constant is obtained from (2.8) by observing 
that since Z^p is independent of p, and = Z7^T^'^\ then 



d 2F(2) 
Up2L g 



p2 

^=ai,T=a2 



= z: 



(2.12) 
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Similarly (2.9) and (2.5) imply that 



Z ,n = ZjT 



(2,1) 
B 



SP=ai,T=a2 



(2.13) 



whilst (2.10) and (2.5) imply that 



Z^=Z. 



.(4) 
2^-3 



As 



(2.14) 



SP=ai,T=a2 



The RG equation can be viewed as a simple consequence of the fact that the bare 
theory is independent of the arbitrary renormalization scale n at which we choose to define 
our parameters. Thus 



(2.15) 
Using the relation between the bare and renormalized vertex functions (2.15) becomes 



K- 



dn 



where 



-lct> 



Za, dK 







Equation (2.16) has the formal solution 

rW(«) = e-^^«>^^"^^rW(«o) 



(2.16) 
(2.17) 

(2.18) 



We can further expand (2.16) by noting that there is implicit dependence on n through the 
renormalized couplings, and possible explicit dependence also. The differential equation 
then becomes 



dn dX 



A2t— - -y 



with 



and 



where 



dt 



.2 = 



N + 



dZ 



-1 



K 



/3(A) = K 



dn 

dX 
dK 



/5(A) 



7A 



IJZ, 



(2.19) 



(2.20) 



(2.21) 



(2.22) 



A Z\ dn 

A being the dimensionful coupling constant. The functions 7^, 7^2 and 7;^ are the Wilson 
functions. Note that a term ^ cannot appear due to the assumed non-renormalization of 
g. If we had chosen to parametrize the anisotropy by a coupling other than a non-linear 
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scaling field, however, (2.19) would then be augmented by a term IgQ^- Equation (2.19) 
in conjunction with the dimensional analysis equation 

P^^=^(^-l)r(") (2.23) 

can be used to relate the correlation functions at two different scales. 

§ 3. RENORMALIZATION GROUP AND THE "ENVIRONMENT" 



§§ 3.1 Comparison of Field Theoretic and "Wilsonian" RG's 

Before proceeding further let us consider the meaning of the RG equation (2.19). The RG, 
in the field theoretic context [31], expresses the invariance of physical quantities under a 
change of parameter. More generally it can be viewed as a reparamctrization invariance 
of a differential equation and its boundary conditions [32]. It is an exact symmetry and 
therefore is a true invariance group of field theory. It is not unlike other reparamctrization 
invariances. Viewed passively it expresses a triviality. Its power is that it can be used in 
conjunction with (2.23) to express a relationship between a parametrized system at two 
different scales. In contrast the RG transformation developed by Wilson and Kadanoff 
[33] is a mapping between probability distributions [34], realized as a mapping between 
renormalized Hamiltonians, and implemented by a "coarse graining" procedure such as 
momentum shell integration [35]. This Wilsonian RG is in fact a semi-group and only 
asymptotically represents an invariance of the system. That the two different approaches 
have led to essentially the same conclusions is a highly non-trivial result. The advantage 
of the Wilsonian methodology is that it is conceptually very clear, whilst the advantage 
of the field theoretic approach is that it is much easier to implement systematically, and 
perturbatively. The general motivation for using an RG approach is as a tool in under- 
standing how systems with many degrees of freedom behave differently at different scales, 
one of the most fundamental facts about our universe being that it is full of "scales" . 

In the Wilson approach it is well known that there are many different types of RG, 
e.g. momentum shell integration, block spinning, majority rule etc. (see [36] for a recent 
review of the Wilsonian approach.) They are all realizable on some sufficiently large space 
of probability distributions, and correspond to maps from measures to measures. Usually, 
however, the Wilsonian RG is taken to be a map from Hamiltonians to Hamiltonians, and 
realizable as a fiow on a sufficiently large parameter space M.. The general intuitive idea 
behind these transformations is that they are mappings (almost invariably approximate) 
between different effective degrees of freedom, represented by an effective Hamiltonian as- 
sociated with different scales. If, for a particular system, each different scheme could be 
implemented, then universal quantities ought to be independent of the particular scheme 
used. An example of such a quantity would be a critical exponent, which is associated 
physically with a system that is almost scale invariant and is determined by a point in M. 
which is associated with the scale invariant system (usually scale invariance also implies 
conformal invariance [37]). If the action of the RG operator is to be a "good" representa- 
tion of the action of the dilatation operator then the fixed points of the RG transformation 
should be coincident in M. with the points of scale invariance. That is not to say that 
they must have the same coordinates, but that the intrinsic geometry of the RG flows and 
the dilatation flows should be the same. For instance, if there were five points of scale 
invariance in whereas the RG used only had three fixed points, then one could not 
think of the latter as being a very successful representation of scale changes. Additionally, 
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the properties of the fixed points should be independent of the "coarse graining" procedure 
one uses, i.e. of the particular choice of RG. 

Many systems, as a function of the parameters of the Hamiltonian describing the 
system, can exhibit scale invariance at more than one point of Ai. Generally one then 
expects to see a crossover between different asymptotic regimes of the theory as governed 
by the various fixed points, hence, one is generically looking at the interpolation between 
two or more scale invariant field theories. If one considered a "coarse graining" procedure, 
such as momentum shell integration, for a system which possesses another length scale 
g~^, other than the correlation length ^, it is clear that the qualitative nature of the 
iterations of the RG should change as one considers momentum shells with k <^ g and 
k ^ g. If one starts iterating at a scale A ^ (7, then as the iteration proceeds into the IR, 
the RG flow will pass close to the g — fixed point (it will only actually hit it for g = 0) 
before proceeding on to the g = 00 fixed point. What is manifest is that a "physically 
sensible coarse graining" procedure will show up the qualitative change in the effective 
degrees of freedom as one considers different regimes of the iteration. The action of the 
"physically sensible coarse graining" follows in a reasonably faithful fashion the action of 
scale changes, hence one more generally may reasonably expect that there exists an RG 
that can successfully describe the change in behaviour of systems as one changes scale, mass 
scale, temperature scale, momentum scale etc. That the Wilson RG has been immensely 
successful is beyond question. Our purpose in the following is to construct such a field 
theoretic RG. 

§§ 3.2 Massive Versus Massless Renormalization 

Let us now turn back to the question of the field theoretic RG. If one considered a massive 
A(/)^ theory in three dimensions, one sets up the field theoretic RG formally in the same 
manner as shown in section 2. If one wishes to investigate the theory under changes in 
scale, and capture the points of scale invariance in then one should look for fixed 
points of a suitable RG. There immediately arises the important question of what set of 
parameters will one's RG transformation depend on. This was not such an issue in the 
Wilsonian approach. There, in principle the RG transformation naturally depends on all 
the parameters of the effective Hamiltonian, in principle an infinite number, though usually 
it is sufficient to only consider it as a function of the relevant and marginal couplings, and 
these are usually finite in number. 

In perturbative renormalization in field theory the story is somewhat different. In 
this setting there is a division of the role played by the Wilsonian RG transformation 
into two parts, represented by equations (2.19) and (2.23). If one accepts the philosophy 
that renormalization is solely a way of accounting for UV divergences, then one has a 
great deal of freedom as to what parameters one's counterterms should depend on, as the 
counterterms are effectively released from any dependence on IR scales. The idea then is 
to find the simplest form of (2.19) by making use of this freedom. The ultimate version 
of this type of approach is minimal subtraction, where one chooses as counterterms only 
those parts of diagrams that survive in the extreme UV limit. For the Hamiltonian (2.1), 
for instance, one knows that in the extreme UV counterterms that are mass independent 
are as good as mass dependent counterterms, if one's purpose is solely to perturbatively 
take the UV cutoff A to infinity. This is because mass is an irrelevant coupling in the UV. 

In our three dimensional example, only two fixed point are perturbatively accessible 
using mass independent schemes, the Gaussian and Wilson-Fisher fixed points with zero 
mass. If one uses mass dependent schemes, another fixed point becomes accessible — the 
infinite mass (temperature) Gaussian fixed point [38] . In fact these statements go beyond 
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perturbation theory — if mass independent renormalization is used then the infinite mass 
Gaussian fixed point cannot be seen from the RG fiow of the Wilson functions, if mass 
dependent renormahzation is used then it can. In other words, in the mass independent 
renormahzation the infinite mass Gaussian fixed point is not captured by the differential 
generator of (2.19), but relegated to a perturbative analysis of the vertex functions. Of 
course one can do such straight perturbation theory near the infinite mass limit, and in 
that sense there was no necessity to implement the RG. The reason one can get away with 
it, in the case of a mass operator, is that the critical exponents associated with the infinite 
mass Gaussian fixed point are mean field exponents, and therefore there was nothing extra 
to exponentiate, however, one is not always so lucky. Typically, if one tries to track the 
theory back into the IR having used a minimal subtraction scheme, in the presence of an 
additional mass scale, gf, perturbation theory will break down, and some other procedure 
will be necessary. It is only the g dependent RG which is capable of giving globally valid 
perturbative information. 

More generally, if one has a field theory parametrized by a set of parameters, P = {g^}, 
corresponding to a point in A^, it might occur that different subsets of the parameters, 
relevant for describing the theory at different scales, are taken into one another by the RG 
flow on M.. If one's renormalization depends only on a subset of the parameters, one is 
restricting one's fiow to take place only in a subspace T oi M.. The resultant RG, RGq-, 
depends only on a subset K of the parameters, and the RG fiows take place only on T. If 
any of the P — K parameters are relevant in the RG sense, then the true RG flows of the 
theory, RGj^, thought of as true scale changes, will wish to flow off T into Ai. However, 
the use of RGj- does not allow for such flows. Such a state of affairs would be shown up by 
the perturbative unreliability of the results based on RGq-. If none of the parameters K 
are relevant then there should be no problem. However, one can only say what parameters 
are relevant when one knows the full fixed point structure of the theory! In principle it is 
obviously better to work with RGj^. If a certain parameter was important then one has 
made sure that its effects are treated properly, and if it wasn't then that will come out of 
the analysis. There can be no danger, except for extra work, from keeping a parameter in, 
but there can be severe problems if it is left out. 

The message then is that in the field theoretic context one's choice of renormalization 
can be quite crucial, some points of scale invariance in M being inaccessible with respect 
to certain RG's. One might be rather worried by this, given that physics should be renor- 
malization scheme independent. One must be careful to make a distinction between the 
points in A^, where the system under consideration is scale invariant, and those points 
which are fixed points of a particular RG. They arc not necessarily the same. If one has 
chosen an RG which is a good representation of scale changes then they will be. We are 
not saying that there are conformally invariant systems that can only be accessed utilising 
certain RG's. What we are saying is that there exist RG's which are sufficient alone to 
describe the system. On the other hand, there are others that must be supplemented by 
extra non-perturbative information. For our massive field theory, the Wilson functions 
derived using mass dependent renormalization are sufficient on their own to access the 
mean field fixed point, whereas for mass independent schemes one has to supplement the 
Wilson functions with extra information from some other source, e.g. summing up possibly 
infinite sets of Feynman diagrams, in order to access the fixed point. 

§§ 3.3 The RG Should be Environment Dependent 

We can now perhaps take a somewhat broader view. In the above we discussed the fact 
that there were inequivalent field theoretic RG's, in the sense that the sets of fixed points 
of these RG's were not the same. Intuitively why is this, and why, in particular, does it 
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happen in the field theoretic RG? Consider a physical system as we observe it at different 
scales. If we start at a very small scale A~^, characterizing the system by a set of degrees 
of freedom (bare parameters), then try to describe the effective physics at some much 
larger scale k~^, it is almost invariably true, in a system with many degrees of freedom, 
that the physics at the scale k is very complicated in terms of the physics at the scale 
A. More often than not a better prescription is in terms of effective degrees of freedom, 
more appropriate for scales ~ n. This is the whole raison d'etre behind effective field 
theory. Of course, experimentally one will find that the effective degrees of freedom are 
very different at different scales. In the above, the effective degrees of freedom at scales 
K ':$> m are essentially massless with power law correlations, whereas for <^ m they 
are essentially decoupled. This is a statement about what one would see experimentally, 
different sorts of correlations in different asymptotic regimes. Clearly one would like one's 
calculational schemes to capture this experimental fact — different effective degrees of 
freedom at different scales. Though the Wilson RG was originally developed to tackle 
systems with just the opposite problem — effective degrees of freedom that are the same 
at all scales (!) — a "good coarse graining" procedure will pick up the change in effective 
degrees of freedom. Similarly, we will show that a "good" field theoretic RG can follow 
the changing effective degrees of freedom. 

We mentioned above that the field theoretic formalism is calculationally much simpler, 
but have gone into some detail about potential pitfalls in using it. In the rest of the paper 
we will see these pitfalls discussed in a much less heuristic fashion than here and also 
see how they can be avoided. The goal then is clear: to develop field theoretic RGs who's 
fixed points encompass all the points of scale invariance for a particular system. The above 
discussion of massless versus massive field theoretic RGs is just an illustrative example, our 
main concern in this paper will be dimensional crossover, however, similar considerations 
hold much more generally as we will now discuss. 

What one considers to be fiuctuations in a system are not unique. Fluctuations are 
usually defined with respect to some "background", or in the terminology we will adopt 
here, some "environment". The environment can have many different properties. One of 
the first things one might ask is what space (or spacetime) does the environment live in. 
Here we want to be more general and consider the space itself as part of the environment. 
Thus if one works in infinite flat space in three dimensions, R"^, or in a "box" of size L 
{L^), then the properties of the space should be counted as part of the environment. The 
reason is obvious. Fluctuations of a field are qualitatively different in and L^. If one's 
only concern is UV divergences then certainly fluctuations with momenta k 3> are 
qualitatively the same on either R"^ or L"^, however, for momenta k ~ the qualitative 
difference is great, and for k -C it is absolutely profound as such modes don't even 
exist on Clearly then the effects of the environment are very scale dependent. If one 
renormalizes one's held theory on using only minimal subtraction say, then, the infor- 
mation one gets out from the consequent RG will be information appropriate to the fleld 
theory on R"^, not L"^. Hence we will say that minimal subtraction is not environmentally 
friendly on L"^, the environmental variable here being L. One could have, of course, an 
anisotropic box with different size sides, then one would have three environmental vari- 
ables. A proper description of the fluctuations from a RG standpoint would demand that 
the RG depended on these three variables. We could consider a space with curvature (con- 
stant or not), this would add additional environment dependence. This would be of some 
relevance in the early universe, where the cosmological constant can also be regarded as 
an environmental variable. The IR behaviour of quantum fields in infiation would only be 
accessible using RG methods using an RG that depended on the cosmological constant. 
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So, one of the most important types of environmental variable comes from properties 
of the space one works in. A second type comes from what sort of background fields are 
present. A gravitational background field is obviously related to the above question of 
what space one is working in. One could also consider the fluctuations of a charged field 
say, coupled to a background electric or magnetic field. The spectrum of fiuctuations in a 
uniform magnetic field B, for instance, is completely different to that for B — 0. The B 
field sets a length scale, and fluctuations with wavelengths very large or small relative to 
it are very different. In fact one can think of the B field as making the fiuctuations live in 
a "box" [39]. The particular background could of course be much more complicated than 
a simple constant field. One example would be the Abrikosov flux lattice, a background of 
monopoles another, the latter possibly being of great relevance for QCD. In other words 
when one is looking at QCD, then the way quarks and gluons behave in a background 
of monopoles will be very different to their behaviour in no background. The standard 
renormalization of QCD, using minimal subtraction, is really only appropriate for the 
situation where there is no background, or one is looking at length scales much shorter 
than the characteristic scale of the background, i.e. at high energy — asymptotic freedom 
— where quarks and gluons are a good representation of the effective degrees of freedom. 

In the above we have given some examples of environment and maintained that the 
environment on certain length scales will greatly influence the nature of the effective degrees 
of freedom which are a description of the system. Quite often the environmental variables 
can be thought of as setting a flxed length scale, or set of scales, which lead to fluctuations 
in the presence of this environment that are qualitatively different at different scales. In 
this paper we will be exhibiting a RG methodology which can describe such changes. 
However, the above example of QCD illustrates another important point, and that is that 
in many cases in practise, one might not be able to make a very clean distinction between 
fluctuations and environment as they interact dynamically. Not only does the environment 
influence the fluctuations but also the fluctuations feed back into the environment. Quite 
often this feedback is negligible, however, when it is not then there needs to be a feedback 
mechanism between the RG as a function of the environment and the environment as a 
function of the fluctuations. A simple example of such a state of affairs is A^^ below the 
critical point. In such a situation the fluctuations influence the background field and one 
can investigate this explicitly via the equation of state. It is the equation of state that 
tells one how the fiuctuations are reacting back on the environment, the environment in 
this case being the background magnetization. For more complicated systems there will 
be more "equations of state", these must be solved in conjunction with the environment 
dependent RG equations to get a closed system. Obviously more often than not this will 
not be easy. 

We hope the idea of environmentally friendly renormalization is understandable and 
hope the reader by the end of the paper will see the necessity of it. 

§ 4. FORMAL PERTURBATION THEORY 



§§ 4.1 Determination of Renormalization Constants 

In this section we will present formal diagramatic series for the correlation functions and 
the renormalization constants. We will assume the generic normalization conditions of 
section 2, but for now will not specify the precise parameter/environment dependence of 
them. The following considerations will be largely independent of such things and therefore 
applicable to a large class of crossover problems. The expressions apply to crossovers where 
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the number of symmetry components of the order parameter doesn't change, however, the 
changes needed to encapsulate this larger class of crossovers are small [18] . 

The effective action (free energy) (2.3) is given by 

r = s + E 

where S is the classical action (mean field Hamiltonian) and S represents the contribution 
due to fluctuations. The diagramatic scries for the effective action can be constructed in a 
loop expansion which gives a series ordered by the number of loops. For the 0{N) model 
under consideration the quadratic part of the Hamiltonian is defined by an operator 

A""^ = AiPf + A2Pf (4.1) 
where the two projection operators Pi and P2 are 

pab = ^B^ and Pf = <5«^-^ (4.2) 



We denote the inverse of 



= GiPf + G2Pf 



For the case of dimensional crossover which we will consider in section 11 

A^^ = {U + ml + tB + ^<P'b)P? + (□ + m| + + ^<p'B)Pf (4.3) 

which gives 

Gi = r and G2 = ^ (4.4) 

n + ml + tB + ^4 n + ml + tB + 

where □ is the Laplace operator in the geometry of interest. The effective action to two 
loops is given by 

T^j cZa;^[0B(n +m| + tB)0B(a;) + ^4(0;)] + ^Tr ln[yli] + i(Ar-l)Trln[yl2] 

+ dx[3Gl{x,x)+2{N -1)Gi{x,x)G2{x,x) + {N'^ -1)GI{x,x)] 

-^JdxJ dy4>B{x)[3Gl{x,y) + {N-l)Gi{x,y)Gl{x,yMB{y) 

(4.5) 

The contributions to the respective F^^^ can be obtained by differentiating with respect 

to (psi^) then setting (pB = (pB^ where (pB is the solution of t'^^\(Pb) = Hb- The F*^^'^'-^) 
can further be generated by differentiation with respect to tB{x). Specializing now to 
constant magnetic field Hb', the magnetization ^b will also be constant for a homogenous 
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and isotropic environment, which we wiU now assume for simphcity. The effective potential 
(free energy density) to two foops is then given by 



j: = 1(^1 + tB)fB + ^4 + ltr\n[Ai] + ^{N - l)trln[A2] 



I dx[3Glix,x) + 2{N - l)Gi{x,x)G2ix,x) + (N^ - l)Glix,x)] (4.6) 
/"^^Z dy^B[3Gl{x,y) + {N-l)Gi{x,y)Gl{x,yMB 



where V is the volume of the system and tr = 

We will now restrict our considerations to situations where 0^ = 0, which corresponds 
to being above the critical temperature with = 0, once again this in no way restricts 
the generality of our methodology [29]. In this case we also have Gi = G2 = G, and the 
vertex functions are proportional to 0{N) symmetric tensors. We will concentrate on the 
coefficient of the relevant tensor eliminating the tensorial structure altogether for simplic- 
ity. Thus adopting a suggestive diagramatic notation (see appendix A for an explanation 
and for the actual expressions in the dimensional crossover problem) and working with 
dimensionless diagrams by pulling out powers of an arbitrary scale k, one finds the one 
loop contributions are given by 

ri = ^»c''0 (4.7) 

rf ' = i^x^.^-^o (4.8) 

rW = _^(^)A|.^-^0 (4.9) 



The two loop contribution is 



^(2) _ :\b ..2d-6 
36 

and 



jr2 = ^iV(iV + 2)«2d-4Q 
A? 



r^-^ = __^K^«-o(iv + 2) {{N + 2)00 + 20 ) (4.10) 



r^"*) = ^^2(d-4) ^^^2 + 6^ + 20)0^ + 2{N + 2){N + 8)00 + 4(5iV + 22)^) (4.11) 



Putting the zero, one and two loop contributions together, we are in a position to 
construct the transformation to renormalized parameters, however, our main concern here 
is to generate the renormalization constants so that the Wilson functions may be calculated. 
An important observation is that for an amputated diagram it is only the number of 
insertions that is relevant not the operator inserted. We have not explicitly exhibited 

the different momentum combinations that contribute to T^^'^ as separate diagrams, the 
sum over the different combinations of momenta being implicit. The diagrams as written 
are therefore to be understood as averaged over the Mandelstam variables, which at the 
symmetric point all become equal. Thus the diagramatic series formally take on their 
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symmetric point form. The reader should however be careful to remember this point when 
implementing expressions for the formal series at a non-symmetric point. 

Adding the different loop contributions together gives 

T = (O + ^nd-\N + 2)0 (4.12) 



rg) =p^ + m% + tB + ^{N + 2)«^-2 ( O - {{N + 2)00 + 2o ) ) (4.13) 



d-4 



.(2,1) _ , _ ^^(AT + 2) (O - ^^^-4^ 



= 1 - (^ + 2) ( O - -^n^-\N + 2)00 



(4.14) 



^ = l--fK2('^-4)(iV + 2)0 (4.15) 



and 



dp'^ 18 

= As - ^(iV + 8)n'-^ (O - ^^^-4(iV + 2)00 

^^^2(d-4) ((7v2 + 6A^ + 20)OV4(5iV + 22) 
36 V 



(4.16) 



where we have regrouped the two loop tadpole terms with the corresponding one loop 
terms. 

The Wilson functions are derived from the renormalization constants Z^, Z^i and 
however, as mentioned in section 2 an additive mass renormalization must also be 
performed, which corresponds in critical phenomena language to a shift in "reference" 
temperature from the mean field critical temperature to some other, usually the "true" 
critical temperature, if such exists. We begin with 

rg^(p, m\ + Ib) =p'^ + m% + tB + ^^b^{p, m| + t^) (4.17) 

(2) 

where is the contribution from fluctuations. There are two undetermined parameters 
in r^'^), m/jg and the scale of Ib- We flnd these by imposing two normalization conditions. 
The first to be determined is the additive renormalization given by = nij^ + Srrij^, 

where m? is the value of T^'^^ at t = and p — 0. We will in fact take m? — 0. Therefore 
we determine by requiring that 

(5m| = sg^O, 5m|) (4.18) 

An important point to make here is that both and Sirij^ are p and t independent, a 
fact that we will use to our advantage later. Condition 4.18, perturbatively for d < 4, is 
very difficult to treat to the presence of IR divergences. The requirement that r(2) be zero 
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when t = and p — can only be used, of course, if the system exhibits critical behaviour 
in the specific environmental setting under study. We then have, as usual, that is a 
linear measure of the deviation from the critical temperature. 

Recall that 

rg\p, ml + ts)^Y: ^r(2-^^)(p, ™|)tf (4.19) 
M 

This implies that 

= rg'0)(0,m|) 

which is our specification of 5m^, however, we can equally well write (4.19) in the form 

M 

where we have just chosen to Taylor expand around a different point, A here being some 
fixed scale. This latter equation implies that at = ^ 

rg)(o, ml + A2) = rg'0)(o, ml + A2) 

All this says is that we can choose any fiducial (RG invariant) point to measure deviations 
in temperature/mass with respect to. It just so happens that if there exists a critical 
point the critical temperature is in some sense a preferred point around which to measure 
deviations. 

Another condition which could be used to determine 5m?Q is 

r(2)(0,t(«e) = «2) = «2 (4.20) 

which corresponds to = 0, Q!2 = Here the scale Kc is a fixed scale being the value of 
the renormalization scale which satisfies the equation t{Kc) = k'^, t{Kc) being the running 
temperature evaluated at the scale Kc- Equation (4.20) gives Sml to be a solution of 

Z^'^kI = + 6ml + ^^■^i^l + ^(0' ™^ + <^™B + ^^^f^l) (4-21) 

If there exists a true critical point, then as the critical temperature is approached, one 
finds that, as Kc — > 0, the above condition becomes equivalent to (4.18). 

The condition (4.18) is to be solved perturbatively to yield the "shift" from the mean 
field critical temperature to another critical temperature. The particular temperature one 
shifts to depends on the values of the environmental parameters, such as ^, entering the 
normalization condition. This is natural because generically depends on environmental 
variables, therefore as the environment changes, the relevant Tc to which it is appropriate 
to shift to changes accordingly. The shift, although intersting, is not something we will 
consider in this paper. The upshot of all the above is that OO (4-14) and (4.16) cancel 
with the mass counterterm associated with the propagators in the one loop diagrams O 
and O- III f^ct a similar cancellation will occur order by order in perturbation theory. 
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Using the condition r(^)(p = 0, t = 0) = gives 

y2 (4.22) 

- -f (iV + 2)/^2'^-6 {(AT + 2)(00 - OloO lo) + 2(o - O |o)} 

where the subscript indicates p = tj^ = 0. However, as mentioned previously this has 
potential problems with IR divergences. An analogous expression which is free of this 
problem follows from (4.21). 

r^"* and r^'^'' at the generic normalization points of section 2 are 



.(4) iN + S) ,2.J-4, 

6 



B =^B 7i O 



(4.23) 



9 ^ ^ 36 ^ ^ 

and 

rg.^) = 1 - (^A,.^-0 + ^-^XWi'-^> (O + (^O^) (4.24) 

Remember that the diagrams in these expressions are now to be evaluated at the appro- 
priate normalization point, explicit expressions for the case of a layered geometry can be 
found in the appendices. 

Differentiating (4.22) with respect to p^, then setting the condition (2.8) yields in the 
two loop approximation 



From (4.24) one finds 



(4.26) 



\2 



36 ^ V/ Q B 



and hence 



(4.26) 



The renormalized four point function, including wavefunction renormalization yields 



+ ^,2(^-4) ((5^ + 22)^ + (iV^+W + 20) + (iV + 2)^^ 
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therefore 

It is obviously trivial to obtain the renormalization constants in terms of the renormalized 
coupling. We leave them in terms of the bare coupling as this coupling that must be kept 
constant when evaluating the Wilson functions. 

§§ 4.2 Diagramatic Expansion of the Wilson Functions 

Having derived the renormalization constants we are now in a position to calculate the 
Wilson functions. Prom the definition of 7^,2 (2.20), using (4.26) one finds 



A2 = — 



(iV + 2),^ ^x„rf-4^^ 

(4.29) 



6 

(iV+2) 



■(^(4-d)0)«''-'A 



y d J, _ d 



6 

where we have defined 

dK dK 

and we now work in terms of the renormalized coupling A. Similarly, from the definition 
of 70 (2.17), and using (4.25) one gets 

M = ^^(^2(4-.)* )'''(''-*'a2 (4.30) 

and finally from (2.21), using (4.28) one obtains for 7;^ 

7A=-^(^,4-<i,0)«<^-*'A 

(5A^ + 22) 1 ~ tAT^o-\ \ . (4-^1) 



+ 
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-^2(4-.)(^ - ^O^) + ^^2(4-.)^) -^('^-^^A^ 



From the definition of 7;^ and its relation to the P function one obtains the /? function 
of the dimensionful coupling 

dX _ (iV + 8) r^)M-4).2 

/(57V + 22) I . {N + 2) ^^ 2(d-4),3 ^"^"^^^ 

+ I ^ ^2(4-d)(^ - 2O ) + ^ ^2(4-d)0 ^A 

We now introduce the dimensionless coupling constant A = Xk'^^^^ the (3 function of which 
is 

/(5iV + 22) 1^2. , {N + 2) \ 3 ^""-''^ 

+ ( g ^2(4-d)(^ - 2O ) + g ^2(4-d)^ ) A 
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We stress that these expressions for the Wilson functions are totally finite for d < 4. 

We now introduce another coupling [27] by performing a coordinate transformation 
on Ai. The new coupling h , which we dub the floating coupling, is deflned so as to make 
the coeflicient of the quadratic term in P{h) unity. From (4.33) this is achieved via the 

redefinition h — — A^_^Q^X. Thus the /3-function for this coupling is 

P(h) = -e{0)h+h^ + ^ ^ '-^ ^ . ^ — ^-h^ 4.34 

where 

s(0)=4-d-«^l^^^ (4.35) 

dK 

Note that in the limit iV — > cxd the /3-function for this coupling becomes exact at the 
quadratic order. The condition P{h) = defines what we call the fioating fixed point, 
the significance of which will be seen shortly. Note that in fixed dimension fiat space for 
massless propagators, or where the mass is taken to be the scale of our normalization 

point, ^4_rfO ^^^d ^2(4-d) (i) ~ ^O^^ ^^e just numbers. Finally we introduce for later 
convenience two combinations of diagrams, fi and /2 defined by 

/2 = (4.37) 



§ 5. ENVIRONMENT DEPENDENT RENORMALIZATION 



§§ 5.1 "Running" Mass/Temperature or Momentum 

In sections 2 through 4 we have set up the machinery for investigating crossovers in a 
formal fashion and without specializing to a particular one. At this point we want to put 
some fiesh on the bare bones of the formalism. Remember that our ultimate goal is to be 
able to describe systems that exhibit qualitatively different effective degrees of freedom at 
different scales. This necessarily implies that the effective degrees of freedom are sensitive 
to the environment. In section 3 wc had a preliminary discussion of this in the context 
of comparing and contrasting the field theoretic RG with the Wilson RG. The previous 
sections were the embodiment of the former implemented perturbatively. So how does one 
see the concept of different RG's and enviroment dependence coming out here? Obviously 
it is the fiuctuations of the system that are sensitive to the environment and it is these 
very fiuctuations that we are cavalierly calling effective degrees of freedom. It should be 
clear then that it is the actual functional forms of the Feynman diagrams themselves that 
exhibit this dependence. 

So how do different RG's appear? This depends crucially on one's normalization con- 
ditions. If, as discussed in section 2, there exists an important mass scale g, then physics 
at scales k, <^ g and k ^ g will be very different. We said earlier that if one's goal is 
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to describe systems at different scales then the RG is in principle an ideal tool for the 
investigation. A good RG will be one who's fixed points represent all the points of scale 
invariance of the physical system in Ai. If we chose renormalization conditions which were 
explicitly g independent, i.e. we used the generic normalization conditions of section 2 
with = 0, for example 

T('^\p,t,\g = 0,K) 2 =(ai + a2)«^ (5-1) 

^=ai,r=Q!2 

and analogously for the other normalization conditions, then the Wilson functions would 
be found to be independent of g. This we emphasize is not just a perturbative failing. 
If the renormalization schemes used are g independent then so are the resulting Wilson 
functions and that is that. If we use g dependent renormalization conditions, however, 
then the resulting Wilson functions will be explicitly dependent on ^. 

To be more concrete: consider (7 to be a fixed "physical" scale (i.e. it doesn't get 
renormalized as mentioned in section 1). We will choose the normalization conditions 

T^^\p = 0,t = K^,X,g,K) = (5.2) 

ar(2) 2 - 

-^{p,t = K ,X,g,K)\p=o^ 1 (5.3) 
r(^)(p = 0,t = K^A,£f,K) = AK^-'^ (5.4) 

r(2>l) (p = 0,^ = ^2, A, £f,K) = l (5.5) 

Condition (5.2), together with the multiplicative renormalization of t, implies that t is 
proportional to T—Tc{g), i.e. that one is measuring temperature relative to the g dependent 
critical point (we are assuming that the system can exhibit critical behaviour for any 
value of g). If an analogous condition with g = had been used, temperature would be 
consequently measured relative to 7c(0), the critical temperature of the isotropic system. 
As these conditions are explicitly g dependent, the RG equation (2.19) for an AT-point 
vertex function takes the form 

(4 ^ -Jm + V(A. f )4 - h^CK f ) {n . 4) ) r<») = ,5.6) 

where now the Wilson functions (3, 7^ and 7^2 not only depend on A, the dimensionless 

renormalized 0^ coupling, in the standard fashion, but also on ^. Notice that 5.2-5.5 imply 
that the parameter which is generating the flow here is the renormalized mass/temperature 
parameter. We could equally well of course have used external momentum to generate the 
flow. 

The RG equation can be solved in the standard manner using the method of charac- 
teristics to yield 



(5 7) 

pK p^K^ 2{pKf ^ Pl^ 
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where A(p), and 0(p); the running temperature, couphng and magnetization respec- 
tively, satisfy 

(5.8) 

dp up 
P^=70^(A(p),^)tW (5.9) 

(5.10) 

With a g dependent renormahzation condition such as (5.4) the P function for the dimen- 
sionless couphng wiU genericaUy have the form 

= -(4 - d)A(p) + a2(^)A2(p) + f; an(^)A-(p) (5.11) 
dp up ^-^ Kp 

Thus we have immediately g independent and g dependent RGs. 

One can now ask how are they different? What is the criterion? Well, let's think of 
the fixed points of the groups. How does one find them in the field theoretic formalism? 
One looks for zeros of the (3 functions. The fixed points of (5.8-5.10) are fixed points of 
the g dependent RG. Fixed points of scale transformations are found by looking for fixed 

points of the dimensionless couplings \,t — and the dimensionless field 6' = ^^^) . 
For the g independent RG the equation (5.11) is 



d\{p) 



00 



, =-(4-d)A(p) + a2(0)A2(p) + J]an(0)A» (5.12) 

The a-n are all just numbers so one gets out a fixed point to lowest order A* = (the 

higher orders can be found in the standard fashion [40] [41]). Remember that a fixed point 

of the RG does not just consist of a zero of one of the /? functions, but actually of all of 
them. In the present case this means one must also consider the equation for the running 
temperature/mass and the running magnetization. Using the renormahzation condition 
(5.1) when, g — 0, implies that t{p) = is a fixed point. One can formally solve (5.9) to get 
t{p) = tcxp J 7^2. Now, t{p) = implies t = due to the multiplicative renormalizability 
of t. The g independence of the normalization condition (5.2) implies that t — T — Tc(0). 
Thus, when = 0, the fixed point being found here is the non-trivial fixed point of the 
isotropic system. An unstable Gaussian fixed point at T = 7c(0), and a stable Gaussian 
fixed point at T = oo, are also accessible to this g independent RG. Clearly no anisoropic 
fixed point is accessible. 

Turning now to the g dependent RG one might think naively that to lowest order 
t{p) = (i.e. T = Tc{g)) and A(p) = ^^(f}\ is a fixed point. However, it is hardly a 

fixed point as the last condition defines a function of p and therefore cannot be preserved 
under the RG flow. Let us change variables to the floating coupling h = a2(^)A(p), the 



(3 function of which is 

CX3 ^ f9_\ 

m = -e{-)h{p) + h\p) + r^TS^h^ip) (5-13) 
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where — 4 — d — ^^j^^^ . To find the fixed points of the system one must remember that 

there is really one more (3 function equation. The dimensionless anisotropy ^ actually 
plays the role of a new coupling in the problem, g. The flow equation for it is 

dg 

^ dp ^ 

The fixed points of all three equations are to lowest order (we omit consideration of the 
Gaussian fixed points) 

^ = 0, /i = £(0), t(^ = 0) = (5.14) 

and 

g = oo, h = s{oo), t{g = oo) = (5.15) 

The first one of course is the same fixed point accessible to the g independent RG. The 
second one, however, is totally inaccessible to this RG. If under favourable circumstances 
one has a good idea of what the cfi'cctive degrees of freedom associated with the g = oo 
fixed point are, then one might impose normalization conditions at (7 = cxo. Once again 
the consequent RG would be g independent, only this time the fixed point accessed would 
be (5.15) not (5.14) as in the case of the g = RG. 

The question naturally arises as to the role of the minimally subtracted RG, which is 
loosely related to the field theoretic e expansion. It is used quite ubiquitously in particle 
physics. The answer is that it is the RG associated with the = 0, t = fixed point. It 
may be used to assist in defining the microscopic theory, in a fashion that is independent 
of an UV cutoff. However, in reality this fixed point plays no more preferred role than any 
of the others. Our formulation makes this manifest. 

The RG that is capable of encapsulating both fixed points is the g dependent RG. It 
is of course by no means unique. Just as in the non-crossover case where different renor- 
malization schemes can exhibit the same physics, e.g. minimal subtraction and massless 
normalization conditions, so here there will be different schemes. The key is that they 
will differ only by relatively unimportant reparametrizations of the crossover variables. 
These differences are the field theoretic analog of the differences between block spinning, 
decimation, blockspinning on different lattices etc. in the Wilson type RG. In the space of 
coupling constants, t, h and g, the physical system exhibits many points of scale invariance, 
two of which are (5.14) and (5.15). One of these is accessible to a 51 = RG, one to a 
g = 00 RG, and only for a g dependent RG are both accessible. In this sense the g depen- 
dent RG is global in the space of couplings whereas the others aren't. The reader should 
understand that we are not saying that using a g independent RG precludes accessing the 
other fixed point, it is rather that the g independent RG on its own cannot access it, 
i.e. the points of scale invariance and the fixed points of the g independent RG are not 
the same. If one is able to supplement the g independent RG with extra information then 
perhaps the other fixed point could become accessible. Almost certainly this will require 
the extra information to be non-perturbative with respect to the g = fixed point. 

§§ 5.2 "Running" the Environment 

In the last section we considered RG fiow generated by mass/temperature or momen- 
tum. Here we will consider briefly the concept of RG flows generated by the environment 
itself. Previously the environment was held flxed as far as the RG flow was concerned. One 
could examine what happened when one changed the environment by examining RG flows 
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with respect to this new changed environment, one could not easily however go directly 
from one environment to another by the action of the flow itself. For the generic anisotropy 
g this means that a particular value of g fixes a curve, one then examines RG fiow along 
this curve. Another value of g yields another, distinct curve along which one can also 
examine RG flows. It is possible however to also generate an RG flow in the "direction" of 
the environment itself. This can generically be done by implementing normalization con- 
ditions with respect to some "fiducial" environment. We will assume that the parameter 
describing the environment does not itself renormalize and therefore there is no distinction 
between bare and renormalized g. For example consider the normalization conditions 

r(2)(p=0,t,A,5r = ^o) = * (5.16) 

ar(2) 

-^^(P,^, A,^ = ^o)|p=o = 1 (5.17) 
r(4)(p = o,t,A,^ = ^o) = A (5.18) 

r(2'l)(p = 0,t,A,^ = ^o,«) = l (5.19) 

The fiducial environment, go^ is abitrary therefore the bare vertex functions must be inde- 
pendent of go. Thus 

where the Wilson functions /3, 7^ and 7^2 depend on A and gQ/\t\'^ ■ The absence of a term 
^99'§g accounted for by the non-renormalization of g. More generally such a term would 

be necessary. 

The RG equation can be solved in the standard manner using the method of charac- 
teristics to yield 



rW(p,,t,0,A,^,^o) = (^oP)'^-^^exp [-^ 



N fP . gox dx 

i<p{Hx),—)— 



r(iV)(i^^J(p) ^A(p)i^,A(p),Al 
P9o {pgor 2{pgoy 9oP 



(5.21) 



where t{p), A(p), and (f){p); the running temperature, coupling and magnetization respec- 
tively, satisfy 

,m=p(X{p),i2^) (5.22) 



dp 



P^ = 7^2(A(p),^)t(p) (5.23) 
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It should be clear here that in this RG, as distinct from that of the last section, it is the 
environment itself that is running. One would implement the same RG strategy however, 
i.e. use the arbitrariness in the scale p to choose a value that allows for a perturbative 
evaluation of the right hand side of (5.21). The appropriate condition to set is ^ = g^p. The 
beauty of running the environment is that, by an appropriate choice of initial condition for 
the characteristic equations, one can relate parameters associated with the environment 
g to those associated with the environment fir = 0. Thus one can answer questions about 
the "shift" in a perturbatively controllable manner using RG techniques. It should also be 
clear that we could combine the two RG's and run the environment and mass/temperature 
together. We will return to these important matters in a future publication. 



§ 6. FORMAL SCALING FORMS 



In this section, working still with an abstract anisotropy mass scale g, we will investigate 
the formal consequences of using a g dependent RG. Consider first the solution (5.7) of 
the RG equation (5.6), where we work now with the floating coupling directly 

rW(t,0,/^,^,«) = («p)^+'^-i^'^ 

RG invariance implies that the right hand side of (6.1) is independent of p, therefore we 
are at liberty to choose a value of it that suits us. We will fix it using the condition 

t{pc) = pIk^ (6.2) 

Pc being the particular value of the RG running scale which satisfies this equation. With 
the generic conditions (2.7) and (2.8), when = 0, one finds = ^ being the 

correlation length. We must then ask which correlation length, as this is sensitive once 
again to what normalization conditions are used, g independent conditions would lead to 
the isotropic correlation length. Here we will assume that g dependent renormalization 
is used and therefore the appropriate correlation length is $,gt, where t — T — Tc{g). We 
append a t to it also because of the assumption that = in the normalization conditions, 
i.e we are using the correlation length in the symmetric (disordered) phase. 

We could have used other conditions to determine pc such as 

^P{Pc)=pW (6.3) 



tiPc) + ^PiPc) = py (6.4) 

The former with appropriate normalization conditions for F^"^^ and corresponds to a 

parametrization in terms of ^g^, the correlation length in the ordered phase when T = Tc{g). 
If fit independent renormalization was used then the corresponding correlation length would 
be that in the ordered phase of the isotropic system at T = Tc{0). The condition (6.4) 
would yield a parametrization in terms of the "true" correlation length of the system, either 
above or below Tc{g). From a purely formal standpoint any of the correlation lengths will 
do, they all correspond to non-linear scaling fields, however, as we shall see, a particular 
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correlation length is naturally associated with a particular field theoretic RG, and within 
that RG certain fixed points of the system might be inaccessible, hence the utility of that 
correlation length from an RG standpoint is diminished. The optimum correlation length 
to use should be the true physical correlation length in the system which is associated 
with the RG that depends on all environmental parameters. Naturally, working with this 
correlation length might computationally be quite difficult. One is always at liberty to 
omit a parameter from the renormalization to try to simplify matters, however, by so 
doing one risks making a particular crossover in Ai inaccessible. Gcncrically a breakdown 
in perturbation theory will then take place, indicating that one is trying to access a region 
of M. which is inaccessible to the particular perturbation theory one is implementing. 

An extremely useful way to parametrize a crossver is via the introduction of effective 
critical exponents [14] , which arc natural generalizations of the standard critical exponents 
associated with one particular fixed point. We define the effective critical exponent for the 
correlation length 

^eff = - -JTZil^ (6-5) 



With the condition (6.2) one sees that 



dln\t\ 



1 



v. 



eff 



2-7</>2 

One can also define an effective exponent for the susceptibility 



dln\t\ 



(6.6) 



where x ^ = T^'^\t,(f),h,g,K) is the inverse susceptibility. An effective exponent r/e// is 
defined naturally via 

c/lnr(2) 
Veff = 2-2- 



(6.7) 



t=0,g 



which is equivalent via the condition (2.8) to 'q^ff = ^(j){g/p)- This 7^, however, is exactly 

the same one that appears in (5.10), only its argument is different. For the exponents 6 
and P associated with the ordered phase one can define effective exponents [29] 



(RnH 



^eff 



and 



dliicj) 
dlncp 



t=0,g 



eff 



dln\t\ 



(6.8) 



(6.9) 



the latter being defined on the crossover coexistence curve. Finally, one can introduce an 
effective exponent for the specific heat 



a, 



eff 



dlnr(Q'2) 
d\n\t\ 



(6.10) 
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There is one last concept we need before proceeding further. Consider in fixed dimen- 
sion d, without g dependence, how F^^) varies with temperature near a fixed point 

dlnr(4) , , , , , 

— = i4-,-2n). (6.11) 

We use here the scahng behaviour of the four point couphng to get information about the 
dimensionality d of the system. One can introduce the concept of an effective dimension- 
ality, dgff, through the natural analog of (6.11) 

dinr(4) 



dint 



= (4 - 4// - 2r7e//)i^e// (6.12) 

9 

dgff is in fact simply related to 7;^, d^ff = 4 — 7;^, and can therefore be thought of as 

a measure of how important the leading irrelevant operator (j)^ is. In the vicinity of a 
particular fixed point 7;^ plays very little role, merely governing the corrections to scaling 
about that fixed point. In a small neighbourhood of this fixed point these are negligible. In 
the case of a crossover this will not be true as it is the corrections to scaling that actually 
interpolate from one fixed point to another. This is why 7;^ takes on an important role. 
In fact, for crossovers where the different asymptotic regimes correspond to systems with 
different upper critical dimensions, d^ff, in the absence of transients, can correctly give the 
change in upper critical dimension as one interpolates between the two fixed points. 

From the g dependent RG equations for the various vertex functions one can derive 
relations between these effective exponents: 

7e// = J^e//(2-r7e//) (6-13) 

^e// = ( 'r^l:':" ) (6.14) 



\deff -2 + T],ffy 

/5e// = ^(4//-2 + r7e//) (6.15) 



O^eff = '2-J^effdeff (6.16) 

and combinations thereof. (RG derivations of these relations can be found in [28], [18] 
and [29]. As long as g <^ k and ^ 3> where k is a typical microscopic scale, these 
exponent relations will be universal, below the upper critical dimension. Just as in the 
standard case where the scaling relations imply that there are really only two independent 
critical exponents, so here there are really only two independent effective exponents. The 
difference with the crossovers considered here is that it is necessary to know one more 
function, 7;^ = 4 — dgff, which represents the effects of the leading irrelevant operator. In 
the case of a crossover for which the dimensionality does not change, e.g. crossover in a 
bicritical system, the effects of the leading irrelevant operator can be subsumed into an 
N^ff, which is a measure of the effective number of components of the order parameter 
[18]. 

Utilizing (6.2) in (6.1) yields 

it, 0, h, g, .) =e|-"-'^e-# /^"^ ^^^rW (1, Xig^,m9^,t)e,t, Kgi,,), g^,,, 

(6.17) 
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where using the running equations (5.8) and (5.10) 

K9^9t)4>\g^gt) = A^^exp /'^(7A - 7^) v (6-18) 
We can rewrite (6.17) in terms of the effective exponents 

r(^) = ) ( 02e-/>"(^«//-2+^^//)^,^^,,) (6.19) 



where ^^^^ is a universal function. Similarly, in terms of one finds 



(6.20) 



In the various integrals in (6.19) and (6.20) one can change variables via the exponent 
definitions, exponent laws, and the conditions fixing p. For example using (6.5) and (6.6) 
one finds 

/ ' (4^/ -2 + 77,//)—= / ideff-2 + VeffHff-r (6.21) 



Thus 

Similarly, one can also find that 

^^/^r_L^ NiJ,„ w^efS di' ,,,, r4> i d^' ''eff d4>' 

= J. i^+d-2id+Veff))j^^^iN)^^^-h ^^^g^-h Jth^) (6.23) 

A particularly interesting case is the equation of state which can be written as 

H = expU 5,ff^jg{te ^-ff^,ge ^eff^) (6.24) 
where ^ is a universal function. The crossover coexistence curve is 

^ = exp| T-^^l (6.25) 



1 Peff 4)' 

Thus using the effective exponents one can derive natural non-linear scaling fields (we 
use the term non-linear here to refer to the fact that they are non-linearly related to the 
linear scaling fields associated with the individual fixed points not as in the sense of Wegner 
that they are eigenfunctions of the dilatation operator). The natural non- linear scaling 

fields that appear above are the various correlation lengths, ^gt = exp j^v^ff{x)^^ ^g^. 



30 



etc. and g (these are also non-linear scaling fields in the sense of Wegner). Note that the 
non-linear scaling field ^gt, interpolates between the linear scaling fields t^'^^^') and t^'^(°°^ 
associated with the g = and g = oo fixed points in the limits (7 — * 0, ^gt 00 {g$,gt ~^ 0) 
and g^gt 00 respectively. In the equation of state there is another non-linear scaling field 



fQ Ji 0eff <!>' which interpolates between the two linear scaling fields ^ifp^Q) and in 
the limits ^ 0, ^ 00 {g^ 0) and g^^ 00 respectively. 

We hope the general pattern is clear. The natural non-linear scaling fields for the 
crossover can be found by taking the natural linear scaling fields for one of the fixed 

points, writing it in exponential form, e.g. = e^i ^ , then replacing the exponent with 

the effective one, e.g. — * e^i "^ff . Note the effective exponent cannot be taken outside 
the integral to recover a form t'^'^ff . In this way one also arrives at the concept of effective 
crossover exponents. Standardly a crossover exponent involves the ratio of two eigenvalues 
of the RG operator linearized around the isotropic fixed point. The generalization of this 
involves exponentiating the integral of the difference of two effective exponents. Naturally, 
one must derive what these effective exponents are. This will be the task of most of the 
rest of the paper. The key is that they should be evaluated using an appropriate RG, i.e. 
one that is appropriate to the crossover in question — one that is environmentally friendly! 

Finally, in this section we would like to discuss universality in the context of crossover 
behaviour. We have stated that in the regime where all length scales are much bigger than 

the lattice spacing a, that the scaling functions etc. and the effective exponents 

will all be universal functions. The meaning of universality here is just the standard one. 
For example, two layered Ising models with different lattice structures will exhibit precisely 
the same crossover curves (up to a trivial constant rescaling) as long as the finite size of 
the system L ^ a. However, the crossover curves for a layered Ising model of size L but 
with different boundary conditions for instance will not be the same. Neither can they be 
made the same by any rescaling. 

The environment affects the infrared behaviour of the theory and therefore the uni- 
versality class. Universality in the non-crossover case is generically governed only by the 
dimensionality of the system and the symmetry of the order parameter. The environment 
gives us extra "labels" for delineating different crossover universality classes. The label 
"boundary condition" for instance allows us to classify, say, three dimensional Ising models 
with one finite dimension into different classes. The effective exponents are different for 
different crossover universality classes. However, the effective exponent laws have a much 
wider universality, being valid irrespective of crossover universality class. For example, 
the effective exponents obey scaling laws even when one of the fixed points is associated 
with the upper critical dimension. Logarithmic corrections to scaling are naturally incor- 
porated in the effective exponents, as indeed are power law corrections coming from the 
leading irrelevant operator. Furthermore, in the simplified context of crossover from the 
Wilson-Fisher to Gaussian fixed points it was found [42] that effective exponents obeyed 
"all the thermodynamic scaling relations (e.g., a^ff + 2^eff + 7e// = 2) but not those of 
hyperscaling" . In fact if one generalizes d in this crossover to 4 — 7;\, we find that all scal- 
ing laws including hyperscaling are obeyed by the effective exponents. Clearly the notion 
of df,ff is inappropriate here, however, the effects of the leading irrelevant operator must 
still be accounted for. This is achieved by the incorporation of 7;^ into the scaling laws as 
above. In the context of the general type of crossover discussed in this paper we have seen 
that the scaling laws are always obeyed. 

We are making an issue of these considerations here for the following reason: often it 
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is the case that the asymptotic scahng behaviour at both ends of a crossover is known, 
i.e. the fixed points and critical exponents at either end are known. Given this knowledge 
Riedel and Wegner [14] wrote down Wilson functions based purely on the requirement 
that the function should allow for an interpolation between the two known fixed points. In 
light of our discussion of universality we regard this as a potentially dangerous procedure. 
To understand exactly why consider the Gaussian fixed point and the Wilson-Fisher fixed 
point for a three dimensional Ising model. These fixed points are the extremal fixed 
points for several crossover systems: three dimensional uniaxial dipolar ferromagncts, three 
dimensional Ising model in a transverse magnetic field, four dimensional layered Ising 
model (with possible different boundary conditions) as well as the standard Gaussian 
Wilson-Fisher crossover of a three dimensional Ising model. If one adopts the approach of 
constructing a crossover curve by simply demanding that the two fixed points are captured 
in the RG flow then there is no way of distinguishing all these physically different crossovers. 
They are induced by different environmental variables which enter the theory in different 
ways. Only via an environmentally friendly RG can one distinguish the different behaviours 
of these systems and see whether they fall into different crossover universality classes. 

In the context of a Wilsonian RG, such as momentum shell integration it is quite possi- 
ble to obliterate real distinctions between crossover curves by implementing the momentum 
shell integrations using too rough an approximation, or by over emphasizing the high mo- 
menutum shells. For instance. Hertz [43] when considering the crossover between quantum 
and classical critical behaviour, derived approximate RG equations which were identical in 
form to those for the Gaussian Wilson-Fisher crossover, by ignoring the relevance of the 
gap in Matsubara frequencies relative to the momentum shell being integrated over, which 
difference is significant at the classical end of the crossover. As mentioned previously these 
two crossovers lie in two crossover universality classes. 

The reader might indeed wonder if the concept of environment is so stringent as to make 
every crossover fall into its own universality class (apart from the standard trivial changes 
in lattice structure etc. one can engender). This is not so. For instance, in the above 
cited examples, it is known [44] [43] [8] that a four dimensional layered Ising model with 
periodic boundary conditions, and a three dimensional Ising model in transverse magnetic 
field, are equivalent. In the language used here they are in the same crossover universality 
class. We will show in section 11 that they ehhibit the same effective exponents. 

§ 7. WHAT SHOULD WE PERTURB IN? 

In the last section we once again worked at a formal level, here we would like to discuss 
the explicit implementation of perturbation theory. We have explained that the main idea 
in crossovers is trying to quantitatively describe systems that exhibit different degrees 
of freedom at different scales. If one cannot solve a model exactly one must resort to 
an approximation procedure. Perturbation theory is a ubiquitous one. However, there is 
always the perennial question — perturbation theory in what? In applications of Wilsonian 
RGs this is often a very difficult question. In the field theory approach it looks much 
simpler, however, our discussion of the different field theoretic RGs should give one pause 
for thought. If one uses g independent renormalization to eliminate UV problems then 
one discovers two things: firstly, that direct perturbation theory fails badly in the regime 
^ 1. New "divergences" appear, which at a given order in perturbation theory are 
generically of the form (gC)^ a,s g^ —>■ oo. Secondly, if one uses a g independent RG one 
finds a running coupling which, in the same regime, becomes very strong. For example, 
for A^^ theory on x one finds 

X{Lk) = 

1 + ^InLn 
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where k is the momentum/mass scale of interest. Clearly in the IR regime, Lk 0, one 
is entering a strongly coupled regime reminiscent of what occurs in QCD. The words "di- 
vergences" and "strong coupling" lead invariably to the invocation of "non-perturbative" 
techniques. It sometimes seems that what is meant by the latter is something which cannot 
be formulated perturbatively in terms of any coupling, such as doing a lattice simulation. 
We believe that quite often the case is overstated. To us, at least in many cases, "di- 
vergences" and "strong coupling" are symptomatic of the fact that one is implementing 
a perturbation theory which is not capturing the correct qualitative nature of the effec- 
tive degrees of freedom. One of the main purposes of this paper is to show that more 
often than not these deficiences can be overcome perturbatively by finding an appropriate 
environment dependent expansion parameter. 

In the case of a theory with an 0{N) symmetry, if one uses a expansion, one finds 

that if a (7 independent renormalization is used, the expansion breaks down in the 
00 limit, except when N ^ 00 before does. This is the spherical model limit, 
which is exactly solvable. One could also consider summing up sets of Feynman diagrams 
such as is done in summing up the daisy diagrams (Hartree/Fock approximation) in finite 
temperature field theory [45]. Naturally, this is always a tricky proposition as certain 
diagrams that are dominant in one regime might not be so dominant in another. One 
must still also address the question of how to renormalize. Additionally, one faces the 
problem that any crossover accessed by the resummation might not be the one of interest. 
This is in fact what happens in finite temperature field theory where the resummation of 
daisy diagrams in an attempt to get an improved description of a phase transition in a 
relativistic field theory merely accesses a mean field fixed point instead of the fixed point 
associated with the transition. 

The key to understanding whether a particular perturbation parameter is suitable 
or not is to see how it gets dressed by fluctuations. If an effective coupling is small in 
one regime there is no guarantee that it will be small in some other. That is why one 
implements renormalization in the flrst place. One could consider a bare coupling to be 
small consistently if one was interested in physics at scales ~ A. However, when one is 
interested in scales p <^ A one finds that the coupling gets large perturbative dressings, so 
one renormalizes to a point k ^ p. In this case one has taken all the fluctuations between 
K and A, absorbed them into A^, and called it A^, which is presumed to be a perfectly 
reasonable, finite, "observable" coupling. The best thing to do of course is consider an 
infinitessimal dressing — this yields the (3 function, the fixed points of which correspond 
to those points where further dressing does not change the coupling. One of the main 
keys then to getting a reasonable perturbative parameter is having a good handle on the 
dressing of the coupling. 

Let us now ask by what fluctuations the coupling is getting dressed? The fluctuations 
that represent the effective degrees of freedom at the scale of interest of course. Usually 
we do not have an exact representation of these fluctuations so we represent them pertur- 
batively. Now, one might think it quite perverse to dress a renormalized parameter with 
fluctuations that were not a good representation of the effective degrees of freedom of the 
system, however, this is exactly what is commonly done, as if for example one uses a g 
independent RG. Such a procedure is not terrible at all scales of course. When k ^ g 
the effective degrees of freedom are effectively g independent and therefore for consider- 
ing physics at these scales it would not be a bad approximation. The fluctuations are 
environment dependent and therefore the dressing of the coupling should be environment 
dependent. If one implements an environment dependent renormalization then one's effec- 
tive expansion parameter is the running coupling which, for example, is a solution of (5.8) 
in the g dependent case. The fixed points of (5.14) are (5.15), neither of these is a good 



33 



expansion parameter throughout the crossover. 

One is not guaranteed of course that the running couphng will be small throughout 
the crossover. This will in fact be the case in this paper where the coupling ~ 1 in certain 
asymptotic regimes. One could take the philosophy that only the lowest order terms should 
be believed in giving the solution, then the higher order terms are used merely to generate 
an iterated form of the lowest order solution. The preferable thing though is to use a 
resummation procedure. For our two loop results we use a [2,1] Fade approximant [40] to 
find 

/3(/i) = -£(0)^ + Y:^ (7.1) 

where 

_ J. ((5iV + 22)A2(4_rf) - ^O') + {N + 2)A2(4_rf)0) 

^" ((A^ + 8)A4_,0)' 

The explicit coupling obviously depends on the crossover in question. Naturally one would 
prefer to work to more than two loops, there is nothing in principle in our formalism besides 
tedious computation to prevent this from being done. So the idea then of implementing 
perturbation theory is the following: one first of all implements it in terms of the running 
coupling constant h generated from a g dependent RG. The actual h used is then the 
resummed solution of the (3 function equation. It is important to realize that in the case 
of the crossover this will be a function, not a pure number, though for a given p and 
which can be arbitrary, it is just a particular number. 

So, what other subtleties does the RG engender? Consider the solution of the RG 
equation (2.16) for F^^) in the disordered phase 

F(2)(t, /^,^, n) = {pnfe- ^^^^^"^ri'^i^, h{p), A, 1) (7.2) 

p^K^ PK, 

Just because one can derive a coupling constant h{p) that is not too badly behaved through 

the crossover does not imply that a perturbative expansion of f(2) will be sensible, the 
reason being that h{p) is only one component in the expansion. h{p) multiplies functions of 

^ and ^2^- Even if h is not too large these functions may become large in certain physical 
regimes, e.g. ^ 0. If this is the case (and this will be the case in crossovers) then one 

must invoke the full power of the RG by using the arbitrariness in the renormalization 
scale p to find a point pc where the functions are not large. If at this point one also 
finds that /i(pc) is not too large then one can make real progress. Obviously pc will be 
explicitly g dependent. The value of pc chosen will usually be determined by the fact that 
one is interested in a regime where the correlation length is large relative to g~^ say. In 
such a regime perturbation theory breaks down in a fairly generic fashion. The thing to 
do, as is well known, is to map, using the freedom of choice in the renormalization scale, 
to a correlation length that is well out of the critical regime. What correlation length 
must one map? We discussed this matter somewhat in the last section. If one uses a g 
independent RG then the natural correlation length is ^q, i.e. the correlation length in the 
isotropic system. Using such an RG, however, will make certain regions of inaccessible. 
The better correlation length to use is the physical one, naturally associated with the g 
dependent RG, which is global and can access all the fixed points of Ai. A good mapping 
point, for instance, would be determined by t{pc) = p^i^^ ^ whereupon pc satisfies 

_..frT<*2(a^) 

.2 



pl = \eK'^A-n (7.3) 
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where t = (T — Tc{g)), i.e a g dependent renormalization has been used. It is this key 
equation which must be solved perturbatively. 

What do we mean by solving it perturbatively? Ultimately what one is expanding 
perturbatively are the characteristic functions of the RG, in particular the quantities 7(^2, 
7^ and 7;^. One determines them as series in the solution of the resummed /? function 
differential equation. One should in fact, based on a democratic treatment of all the 
characteristic functions, Padc rcsum the series for 7^2 and 7^ as well. To two loops, 

obviously 7^ cannot be resummed. The [2,1] Fade resummed expression for 7j,2 is 



A2 = 



^ + 2^ ' (7.4) 



N + 8J e A,(4-.)(^-|01 I 1 -42(4-.)^ ^, 

For small values of N one finds that the effective exponents derived from [2,1] Fade re- 
summed Wilson functions differ little from those derived from a using resummed /^-function 
for the floating coupling with the perturbative series for the effective exponents in this cou- 
pling. See Fig. 10 for an illustration of this in the context of dimensional crossover. The 
former prescription caprures the exactly solvable n — 00 limit and therefore seems prefered. 
Correlation functions are related to exponentials of the Wilson functions, it would be quite 
against the spirit of the RG to start perturbatively expanding the exponentials, otherwise 
there would be little point in using the RG in the first place. This is one of the problems 
with trying to use e expansion methods to calculate scaling functions. The question is 
always: if something is to a power of e should one e expand it? In (7.3) one should expand 
7^2 in the integrand to a given order k"', h being the resummed solution of the P function 

differential equation to order h'"''^^. One should not expand the exponential. To 

lowest order for instance, one would solve 



2 t 
Pc = ^exp 



N + 2\ fP^ , ^ g,dx 
h{x,-) 



N + 8J Ji ^ x' X 



(7.5) 



where /i is a solution of 



p^ = -e{0)h + h^ (7.6) 
dp 

More often than not one will have to resort to numerical techniques to solve the equation 
that fixes p, however, we will see shortly that as long as one captures the essential physics 
of the crossover with an appropriate RG then one's final answers are rather robust under 
changes in solving the condition. 

Besides the ambiguity in using the e expansion to determine scaling functions, i.e. 
should or shouldn't one expand something which is to a power of e, there are two other 
defects. First of all, if one makes a definite choice in deciding to expand everything in £, 
as is done in finding the equation of state [46] , one will find that the resulting expression 
is only thermodynamically valid in a restricted region of the thermodynamic state space. 
Secondly, if one is considering a crossover wherein the upper critical dimension changes 
from one asymptotic region to another, then an e expansion will be inadequate, being able 
to capture the singularity at one end but totally failing to capture the singularity at the 
other. In principle, in such circumstances a expansion could work, if used in conjunction 
with an appropriate environmentally friendly renormalization. However, it is well known 
that the expansion gives poor results for small values of A?", is difficult to implement at 
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higher orders, and additionally such an expansion would be useless if there was a change in 
the number of symmetry components of the order parameter during the crossover, e.g. for 
the bicritical crossover. We believe our perturbation technique to be the only one on offer 
applicable to crossovers where both the dimensionality and/or the number of components 
of the order parameter can change. 

In the above we are exploiting the arbitrariness in the RG scale to allow us to map to 
a point where perturbation theory can be trusted. We now discuss how one can exploit 
the various renormalization conditions in order to facilitate calculations. We have already 
seen that the type of condition used, in the sense of dependence or independence of the 
environment, is quite crucial to achieving a description of crossover systems. There are, as 
pointed out, different, inequivalent representations of the RG. Having determined the RG 
that is most suitable, and having determined the RG scale that allows perturbation theory 
to be implemented across the crossover, one can make further simplifications by trying to 
put as much of the physics as possible into the characteristic functions alone. 

What this means is that for an object like the inverse susceptibility, choosing a renor- 
malization condition like (5.2), implies 

r(2)(i,/i(p,),^,i) = i (7.7) 

which gives 

r(%A,a,.)=exp/'f|^)- (7.8) 

Clearly through the choice of normalization condition all the physics has been put into 

the two functions 7^ and 7^2; i.e. r(^)(l, h{pc), 1) did not need to be perturbatively 

expanded. According to the framework outlined above, a perturbative treatment of (7.8) 
entails expanding /3, 7^ and 7^/,2 as a power series in /?. to a certain order, Pade resumming 

the resulting series, and then putting in h as the resummed solution of the P function. We 
could have used a condition other than (7.7), such as 

r(2)(o,/,(p),A,i) = o (7.9) 

. In such circumstances r(^)(l, h{pc), 1) would also be perturbatively expanded in h, 

hence not all the physics would be in the characteristic functions. One would obtain an 
expression of the form 

T(''\t,h,g,K) = e^p[j:[^)f) (^l + h{pc)B,{-^J+hHpc)B2{^) + --^ (7.10) 

where Bi and B2 are dependent on the precise normalization condition one uses. For 
the massless normalization condition (7.9) Bi contains a contribution from the difference 
between the massive and massless tadpole. At the one loop level we could have exploited 

the arbitrariness in the renormalization scale p to set T^^^{l,h{pc).,-^,1) = 1. In this 
case 702 will not be the same as found using the condition (7.7). A quantity such as an 
effective exponent, however, will be invariant under such a change in normalization point. 
Unfortunately this would not work beyond the one loop level, for d < 4, due to the presence 
of IR divergences in B2, as was pointed out in section 4. If one wanted to calculate a scaling 
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function, such as the susceptibihty above, the functions Bi and B2 would, for instance, 
carry information about the difference between the GMS method used in [15] and [19], 
and using a normahzation condition, as we have used in our method. Although they might 
not be IR divergent, as emphasized above, this difference is not easily exponentiated, so 
that in the calculation of a scaling function GMS will yield less information than using 
normalization conditions. This is one disadvantage of GMS versus our method. 

It is useful to give an analog expression diagramatically to lowest order in the coupling 

A 

r(2)(t, h, g, k) = exp (/;(i_(^(^(,_,,0)x)*) (^1 + ^(AT + 2)(0 ' - O ) + • • •) (7.11) 

The quantity (O ' ~ O ) is a measure of what has been "left out" of the RG and would be 
the analog of ln( Jj) in four dimensions, which is the difference to lowest order between the 

normalization condition T^'^'^ (t) = t and minimal subtraction. The RG equation tells us how 
to perform an infinitcssimal dressing of a quantity. By different choices of counterterm one 
can dress the quantity with different fluctuations, g dependent, g independent etc. Among 
the g dependent schemes themselves though, as we are seeing here, one can dress things 
differently. Should one dress things using the whole of a Feynman diagram, or what we 
think to be the most important part of it, e.g. the mass independent part, or....? The RG 
only tells us about dressings, but because we can integrate an infinitessimal dressing, one 
can, in principle, derive a huge amount of information from it, as this integrated dressing 
gets exponentiated in the process of solving the RG equation. What is not included 
as dressing cannot be exponentiated, therefore one potentially loses a large amount of 
predictive power. Manifestly the best thing do do then is dress up as much as possible, 
and particularly in an environmentally friendly manner. 

§ 8. THE FLOATING FIXED POINT 

In this section we will introduce the concept of the floating fixed point. We said in section 6 
that to find the fixed points of the RG one had to find the fixed points of all the (3 functions, 
including j3{g). We exhibited the two relevant fixed points. If we had just examined /5(A) 

in isolation, and set /?(A) = 0, to lowest order one finds A = '^"^ — -z— . In the limit 

6 (^(4-d)0) 

^ — > 00 one finds that A — > 0. This implies that it wasn't very useful to think of the 

algebraic zero of the (3 function as a fixed point as it did not cross over between the right 
values (a more detailed discussion of this can be found in [27]). If we use the floating 
coupling, however, and consider f3{h) = 0, one finds that the solution of this equation 
interpolates between the correct fixed points of the crossover. The condition f3{h) = 
defines the floating fixed point h* . To lowest order h* = ^(O)- To higher order one could 
try to iterate a solution, i.e take this solution, put it into the h"^ term in (3{h), and resolve 
the algebraic equation /3(/i) = 0. To the next order the iterative solution is 

h*{z) = e{z) + —1-2 ((5iV + 22)h{z) - {N + 2)h) e{zf + 0{e{zf) 

Finally one can Fade resum the floating fixed point to obtain 

h*i \ ^^^^ 
h {z) 
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The utility of using the floating flxed point is that to flnd it one need only solve an 
algebraic equation instead of a differential equation. The floating flxed point is guaranteed 
to give the correct asymptotic answers, i.e. the floating flxed point coincides at ^ = 
and g = oo with the isotropic and anisotropic fixed points which are the fixed points of 
the system of differential (3 function equations. Naturally the solution of the differential 
equation and the algebraic equation will differ through the crossover. This difference is 
found (in all cases examined so far) to be relatively small, so the floating flxed point can 
be thought as giving a good approximation to the true crossover h. 

One can go through the entire analysis of sections 6 and 7 expanding everything around 
the floating fixed point in analogy to expansion around a normal fixed point. For example, 
the non-linear scaling fleld expanded around the floating flxed point becomes 

= A{t,g)e^p (8.1) 

where A{t^g) — exp //(z^e//* — ^eff) is a correction to scaling factor around the floating 
flxed point and v* is the effective correlation length exponent evaluated at the floating 
flxed point. Remember that i^eff* is a function of t and (7, therefore one cannot pull it out 
of the integral as if it were a constant. This can only be done asymptotically where v* 
becomes the exponent at the true flxed points. If we had expanded around the isotropic 
flxed point we would flnd 

ig, = B{t,g)e^p{- (8.2) 

where B{t,g) = exp Ji{i^{0) — i^eff) is a correction to scaling factor around the isotropic 
fixed point. The key difference is that A is small throughout the crossover and approaches 
zero in the asymptotic regimes, whereas B goes to zero near the isotropic flxed point but 
becomes singular near the anisotropic flxed point due to the fact that it must contain a 

singularity of the form in order that the isotropic singularity be replaced by 

the anisotropic one. Thus using the floating flxed point corrections to scaling can be made 
small throughout the crossover. Importantly, in our formalism the corrections to scaling 
around the floating flxed point can also be calculated explicitly. One can deflne effective 
exponents, cte//*, A//*: 7e//*: ^e//*: ^eff* aud an effective dimension d*fj: associated with 
the floating flxed point and flnd that they obey exact analogs of the equations (6.13-6.16). 
Similarly for all the non-linear scaling flelds, there are floating flxed point analogs which 
differ from the one's already deflned only by small corrections to scaling throughout the 
crossover. 

§ 9. DIAGRAMMATIC REPRESENTATION OF EFFECTIVE EXPONENTS 
AND SCALING FIELDS 

In the last sections we tried to convey the utility of the fleld theoretic RG if implemented 
judiciously. In this section we wish to derive perturbative expressions for the effective 
exponents and some of the scaling flelds mentioned in section 6. Things will still be formal 
to the extent that we do not specialize yet to a particular crossover. The reason for doing 
this is that in our diagramatic notation it is clear that any crossover can be treated by 
putting into the diagrams the effects of the particular environmental parameters that are 
inducing that crossover, whatever they may be. We have already given the expression for 
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the Wilson functions diagramatically so we just need to construct the effective exponents 
from them and deff. Note that here we do not Pade resum 7^2 and we implement the 
scaling laws perturbatively in h. These expressions will be useful for comparative purposes 
and as we shall see are quite reliable for small values of N. Analogous expressions can 
be obtained when 7^2 is also resummed, in which case the effective exponents should be 
constructed from these resummed quantities without further expansion. 

One finds to two loops 



Veff 



, (Ar + 2) ^2(4-rf)^ 

\N+8f (A^-dOf 



_ 1 (iV + 2) 
^^^^ - 2 + 4(Ar + 8)^+ 

(iV + 2) 



7e// = H 



8(iV + 8) 
(iV + 2) 



iV + 2 + 12 



+ 4- 



h + 



(iV + 2) 



A 



2{N + 8) 4(7V + 8)2 
Using cZe// = 4 — £(0) ~ one finds 



iV + 2 + 12- 



2(4-d) 



deff = 4:-h- 



A, 



{N+8fy {M-dOf 



(9.1) 



(9.2) 



(9.3) 



(9.4) 



and hence 



(4 -TV) , 
^ft 2(iV + 8) 



f^eff — 



( (4 - iV) (iV + 2) + 16(iV + 8) ^^'^"^^ ^^'^ ^ 



(^4-dO) 



4(iV + 8)^ 



/ 



2(iV + 



3(iV + 2) + 2(7iV + 38) -^^(^-'^)(^-^p) \ 

(^4-dO) 



4(iV + 8)2 



' ^(5Ar + 22)^2(4-d) - lO^) 

1 + 8 7, -F, 



2 I (AT + 8)2 (A4_,0)2 



/ 



(9.5) 



(9.6) 



(9.7) 



If one solves the Pade resummed differential equation (7.1) and inserts it in these expres- 
sions then one has the effective exponents for a general class of crossovers. For the Pade 
resummed floating fixed point 



£(0) 1+ 



■'2(4-d) 



(Af+8) 



^ (5Ar+22) 



-1 



, + (^+2)^^(1^^ I £(0) 



(^4-dO) 



(9.8) 
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Thus one finds 



Veff 



*2(4-d) \V-2 



1 (TV + 2) 



2 ^ + (JV+2) 



^2(4-d)' 



(^4-dO) (^4-dO) 



A. 



2{4-d) I V-2 



(9.9) 

-1 



2 4(7V + 8) 
+ — ^(iV+2+12- 



(jV+8)2 I (5^+22) ^ ^ 2 ^ + (^^+2)- 



(^4-dO) (^4-dO) 



(^4-dO) (^4-dO) 







8{N + 8) 

^ \ (^4-dO) (^4-dO) 



(9.10) 

With the effective dimension at the floating flxed point, d*jj = £(0)i oi^*^ can easily derive 
all the other effective exponents at the floating flxed point. We leave this as an exercise 
for the reader. 

Consider now the non-linear scaling flelds ^gt = e Ji *' and x = te ^ ^^/z <t> _ 
D iagr amatically 



exp 



and 









X = t exp 




l+(jV+8) + 



3(iV+12)+2(7JV+38)- 



^2(4-d) 



(^4-dO) 



2(iV+8)-^ 



(9.12) 



In terms of the Pade resummed floating flxed point 
Cgt = exp 
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1 I iN+2) 



■^2(4-d) I ^~2 
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8(iV+8)2 ^ ^ ^ 



^^+4^^^^).(0)^X 



^2(4-d) 



+(^+2) --'7^2 
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(9.13) 
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and 
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t exp 



i+livT8y^(0)| 1+ 



(N+8y 
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^2(4-d) 
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(9.14) 

It is impressive to see from these expressions just how much diagramatic information 
has been extracted using the RG. The reason we have exhibited these cumbersome looking 
expressions, is that they can serve as "black boxes" for evaluating crossover functions. 
For any crossover falling within the general class treated, in our two loop perturbative 
analysis, the reader may take the expressions given, insert the appropriate propagator 
into the diagrams, crank the handle, either numerically or analytically, and obtain the 
corresponding crossover scaling information. Now it is time to turn to a specific example 
of a crossover to show all the ideas and formalism we have presented at work. 



10. EXPLICIT DIMENSIONAL CROSSOVER RESULTS 



§§ 10.1 Periodic Boundary Conditions: One Loop Results 

The specific crossover we will consider here is that of a d dimensional layered system 
satisfying periodic boundary conditions (or in particle physics language \4>^ theory on 
R^"^ X S^). We consider, to begin with, an 0{N) vector model. We will be concerned 
here with T > Tc{L), where Tc{L) is the critical temperature of the finite size system, and 
d < 4. We will also assume here that the finite system does exhibit a critical point, though 
the formalism is equally applicable to the case where such does not exist. The environment 
here is the thickness of the system, which we denote by L, and the periodic boundary 
conditions. The results below are equally applicable to relativistic finite temperature after 
the trivial substitution L = 1/T and = l/rriT. 

Let us begin by examining the one loop expressions 

P{h) = -e{pKL)h + h'^ + 0{h^) (10.1) 



70 = to one loop. The solution of (10.1) is 



hi^Lp) = — I (10.3) 

41 



Choosing ai = 0, a2 = p% in our normahzation schemes we have pcn = Thus we have 
h in terms of z — L/$^i, and an initial coupling at the scale k, or equivalently zq = kL. 
Substituting this into the one loop contribution for 7^2 we have 

(N + 2) (T L^^^^^^ 

Now, in d < 4 there is a non-trivial fixed point of the bulk system, and hence if we impose 
the initial condition on the RG flow at = 00, i.e. at zero correlation length, and choose 
a finite renormalized coupling h{oo), we will be probing the universal part of the crossover. 
This is not a necessity but merely allows us to isolate this universal part. Often, one is 
also interested in corrections to scaling, if zq is not chosen to be infinite then these are 
automatically included in the expressions, and as emphasized in section 7 the effective 
exponent scaling laws are still valid in this case. One finds 

h{z) = (10.5) 

which depends on only the variable z, and interpolates between the bulk and reduced 
fixed points. This is the equation of the separatrix between the two fixedpoints in the one 
loop approximation. It is not difficult to verify that (10.5) satisfies the one loop (3 function 
equation. Note that the component diagrams are finite for d < 4 but for (i = 4 O diverges, 
it is therefore impossible to eliminate the dependence on zq in this case. 

With periodic boundary conditions we have 

°° 47r^n^ 47r^n^ 



E 



1 + 



2;^ V z"^ 



e{z) = 5 - d - (7 - d)^^=^^5 (10.6) 

^ / 47rV\- 



n=— 00 



^2 



and the universal one loop ffoating coupling 

00 (d-7) 

Kz) = (5 - ^ (10.7) 

E (! + (¥)') " 

n=— 00 

From these expressions the effective exponents can be evaluated. Rather than writing the 
expressions down explicitly, as we will be presenting shortly the full two loop expressions we 
will restrict ourselves here to making some specific comments about particular dimensions. 

Conveniently for d = 3, e{z) and h{z) are expressible in terms of elementary functions. 
We have explicitly 

, , 2;2coth(§) 

e{z) = 1 + ^-r— P 10.8 
smh z + z 
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Similarly an explicit form for h{z) is 

^ zsinH^f I sinhzo + ZQ _ ^cotKf) ^ sinh. _ 
smhz + z h[zQ)zQsinh{^)^ zq smhz + z 

We present h plotted against ln(l/2) in Figure 22, where we choose Hq = .9 for ln(l/2;o) = 
—7.5. The crossover to meanfield theory is evident as ^ — > oo. 

Focusing on the universal part of the crossover by choosing zq —>■ oo, with h{zQ) finite 
we find 

h(z) = l + -^ (10.10) 

smh2; 

Since in the case of the two dimensional end only the = 1 model exhibits a standard 
second order phase transition we restrict our considerations to this model. (The techniques 
used here though can also be used to treat other values of A?^, in particular dimensional 
crossover in a non-linear a model [47]). We thus find 

As 70 = to one loop so is rj^ff- Similarly, substituting for h in (9.2) for v^ff to one loop 
gives 

Asz — > OO, I'eff — 0.58 and as 2; — 0, I'eff — 0.67. FVom (9.3) our expression for 7e// is 

leff = 1 + (10.13) 

6 6smh2; 

which varies between 1.17 and 1.33. To one loop cZg// = 3 — 2;/ sinh2; and varies between 
3 and 2 as 2; varies between 00 and 0. Using the scaling law (6.14) one gets 5e// = 

4+ g^^Yiz ^^i*^^ varies between 4 and 5; and from (6.15) one finds /Jg// — ~ q ginh z ^^^'^^ 
varies between 0.33 and 0.17; and finally a^ff = ^(1 + ^13^)- O^ig can also determine 
approximations to the effective exponents by calculating them with respect to the fioating 
fixed point. We emphasize once again that one can independently verify the scaling laws 
by studying the system when T < Tc{L) (see [29]) and calculating the effective exponents 
directly. 

In the case of (i = 4 we are not as fortunate in finding expressions in terms of elementary 
functions for our one loop results, however, the expressions involving the sums themselves 
are not too bad and are rapidly convergent. Explicitly 

°° 47r^n^ / 47r^n^^ 



, n=— 00 



1 + 



eiz)^!-^'^^^ (10.14) 



1 + 



n=— 00 



The one loop solution to the /3 function is given by (10.3). It is now not possible to set the 
initial value of the coupling to a finite value and take the limit of 2:0 —> 00. The logarithmic 
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corrections to scaling in d = 4 mean that the running coupling will retain a dependence 
on the the initial condition, i.e. the scale k will remain even when L drops out. Since the 
one loop results have been presented elsewhere [27] we will not repeat them here. 

§§ 10.2 Periodic Boundary Conditions: Two Loop Results 

Prom the considerations of the previous sections we can conclude that the Wilson functions 
are given by 

(3{h, z) = -e{z)h + h'^- ^ M hN + 22)h{z) - (iV + 2)h{z))h^ + 0{h^) (10.15) 

(7V + 8)^ 

^'P = 27^r^/2(^)/^' + 0{h^) (10.17) 
(iV + o j 

where the functions /i and /2 applicable for general d for this particular crossover are 
given in Appendix A. In Appendix B the explicit, analytic expressions for the case d = 4, 
«! = 0, 02 = 1 are presented. Recalling that v^ff = 1/(2 — 7^2) and rj^ff = 7^ we have to 
two loops that 

= 1 + + (iV + 2)(iV + 2-12A(.)+4/,(.)) , 

2 4{N + 8) 8(iV + 8)2 ^ ^ 

Similarly, from the scaling relation for ^^ff {'Jeff = ^e//(2— ?7e//)), which we have established 
in section 7, we have 

and from the defining relation de// — ^ — e{z) — we have 

4// = 4 - + ^^((SAT + 22)/i(^) - (AT + 2)/2(^))/i2 (^q ^O) 

\N + 8) 

Using the scaling law for Oe// (cte// = 2 — v^ffdeff) yields 

'•^^ 2(iV + 8) 4(iV + 8)2 ^ ^ 

The scaling law for (3gff {(3eff — -^^^(de// — 2 + rj^jj)) gives us 

. _1 3 (3(Ar + 2)-2(7iV + 38)/i(.)) 2 .^.22^ 

- 2 ~ 2(ivT8) ^ 4(iV + 8)2 ^ 

And finally using the scaling law for S^ff {S^ff = {d^ff + 2 — r]eff)/{deff — 2 + rjeff)) gives 
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Note that the functions e{z)j /l(^) and f2{z) are independent of N. For a three 
dimensional layered geometry {R^ x 5"^), we find e{z) interpolates between 2 and 1, fi{z) 
interpolates between 0.28 and | and f2{z) interpolates between 0.23 and ^ as z varies 
from to oo. For a four dimensional layered geometry (R"^ x S^), e{z) ranges monotonically 
from 1, for z = 0, to for 2; = 00, and takes the value ^ for 2 ~ 3.3. The function fi{z) 
ranges from ^, for 2; = to ^ for ^ = 00; f2{z) ranges from ^ for 2; = to | for 2; = 00. 

Thus the equations (10.15), (10.16), (10.17) interpolate between those for and i?^ as z 
ranges from to 00. We provide plots of s{z), fi{z) and f2{z) in Figure 1 for the layered 

four dimensional geometry with periodic boundary conditions {R^ x S^). 

To explicitly implement the two loop results we resort to a Fade resummation tech- 
nique. As is well known it sometimes occurs that the zero of the /3-function disappears 
at two loops only to return at three loops. This is a property of the perturbative series 
and is not specific to our problem. It is equally true in the fixed dimension non-crossover 
case treated by Parisi. We therefore adopt Parisi's approach and use a [2,1] Fade resum- 
mation. As this agrees well with high temperature series and s expansion results there is 
every reason to have confidence in it as a good method of capturing the true nature of 
the resummcd expressions. For three dimensions the two loop Fade resummed results are 
in excellent agreement with six loop resummed results [48], where the Callan-Symanzik 
equation in distinction to our homogeneous RG was used, and the best high temperature 
series [41] . A Borel transformation is also an option, we have not as yet followed this route. 
Our series are still too short to make this worthwhile. Calculations to higher orders could 
be implemented, and it should not be too difficult in the case of a three dimensional layered 
geometry to adopt the numerical techniques of [48], in conjunction with our diagramatic 
approach, to getting numerical curves for the expressions to a similar accuracy. 

We have also resummed the series for 7^2, in terms of the Fade resummed coupling 
and used the scaling laws for the effective exponents. The differences between these results 
and those obtained when 7^2 is not resummed, and the scaling laws not used, are small 
for small values of N. We plot in Figure 10 a comparison of the one loop, two loop 
Fade resummed 7^2 and non- resummed 7^2 results for ^e// with N = 3. By adopting the 
philosophy that all the characteristic functions should be treated on an equal footing, in 
other words they should all be Fade resummed one retains the spherical model limit as an 
exact limit of the Fade resummed expressions. After Fade resumming of the characteristic 
functions the scaling laws should be used directly without further perturbative expansion. 

Our procedure is therefore to solve the Fade resummed differential equation z^ = P{h, z) 
for h{z) and substitute it into the Fade resummed expressions for 7^2, 7^ and use the 
scaling laws to obtain the effective exponents. We were forced to do this numerically as 
the (3 function equation does not have an analytic solution known to us. 

An alternative procedure would have been to solve the differential equation iteratively. 
The one loop equation is solved first, and the solution substituted into the two loop portion 
of the (3 function equation. The resulting differential equation is then solved. The pro- 
cedure is iterated as successive loops are calculated. We have not adopted this approach, 
however, it would be interesting to compare the accuracy of the two procedures, this could 
only be reliably done at higher order than two loops. 

The condition P{h, z) = defines from section 7 the floating fixed point, valuable 
in describing a system where, in terms of the usual picture of flow lines associated with 
the /3-function, the flow lines are explicitly "time" dependent. Again note that in the 
ai = 0, 0:2 = Pc'^^ case z = pcuL. The floating fixed point, as a solution of an algebraic 
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and not a differential equation, is manifestly independent of zq. In a comparison of the 
running and ffoating effective exponents, as the values of z may differ by a constant scale 
difference, one curve may be shifted along relative to another. There is also of course some 
mismatch between the solution of the differential equation and the floating fixed point in 
the crossover region even after this shift is accounted for. This is because solving the /3 
function differential equation is not the same as solving the associated algebraic equation. 
They do, however, track one another closely (when the solution of the differential equation 
is dominated by a non-trivial fixed point) but not exactly. The fact that the fioating 
fixed point is independent of initial conditions means there is no dependence on transients, 
hence logarithmic corrections to scaling are suppressed, however, these can be recovered 
by examining corrections to scaling around the floating fixed point. 

The zero of p{h, z) from (10.15) is 

h*{z) = e{z) + —1-2 ((5Ar + 22)h{z) - {N + 2)h{z)) e{zf + 0{e{z)f 

(iV + oj 

Substituting this value of h into (10.16) and (10.17) gives 

= I^tI'^"^ + ^ (^^ ^^^^ ^ " " ^^^'^"^^ 

and 

These functions interpolate nicely between the two asymptotic fixed points at z — and 
z = oo, and define finite size scaling functions around which one can develop a correction 
to scaling expansion, in h — h*. The analogous [2,1] Pade resummed expression for the 
floating fixed-point is 

h*{z) - ^^"^^ 



^ ~ ^^^^ + ^^^-^i^^^ ~ + ^^^2(;^)) e{z) 



We leave it to the reader to substitute this into the corresponding expressions for the 
effective exponents. Figure 11 shows a comparison of the Pade resummed floating coupling 
and floating fixed point for d = 4 and periodic boundary conditions. 



§§ 10.3 Pade Resummed Two Loop Results for Some Specific Models 

We will now write down explicit results for some models of particular interest. The two 
loop results we consider explicitly here are for general d. The explicit expressions for fi 
and /2 can be found in Appendix A by setting ai — and a2 — 1- Specific expressions for 
d = 4 are contained in Appendix B. As the expressions are rather unwieldy we choose not 
to exhibit them here. We should also point out here that these two loop Pade resummed 
results hold identically for the corresponding finite temperature relativisitic field theories. 
We present the direct perturbative results for ease of comparison with other methods. 

For the Ising model, AT = 1, we have 

Pih, z) = -e{z)h + h^-^ (9/i(^) - f2iz)) + 0{h^) 
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Ah, (/l(^) - Im^) ) + 0{h') (10.24) 



MKz) = -f2{z)h' + 0{h'') 
We find the resulting effective exponents to be 

_ 1 1 . , (3-12/i(z) + 4/2(z)) 2 

leff-l + ^h+ h 

d^„^4-h+'-^^M^^M^h' (10.25) 

^^■^ 2 6 36 

For A'^ = 0, wfiicli is related to the self avoiding random walk and polymers, we have 
(3{h, z) = -e{z)h + h^-\ (ll/l(^) - f2{z)) + 0{h^) 

Mh, ^) = 1/i - A (^f,{z) - ^/2(^)) + 0{h^) (10.26) 
The corresponding effective exponents are then 

1 ^ (i-oAM+s^^, 

2 16 128 

1 (1-6A(^ 2 

4,/ = 4 - /z + (ll/l(^)-/2(^)) ;,2 (10 27) 

8 

_1 (l-16(/i(.)) 2 

I 3 (38/i(.)-3) 2 
^^^^-2~T6''+ 128 ^ 

... = 3 + .+ ii^.^ 
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For N = —2 which is beheved to have the same two point correlation function as the 
Gaussian model we have 

P{h, z) = -e{z)h + h^- hi{z)h^ + 0{h'^) 

7^2(/i,^) = (10.28) 
l(j,{h, z) =0 

As we see there is a non-trivial crossover of 7;^, this induces a crossover in cte//, /?e// and 
Seff. The other exponents retain their mean field values, z^e// = 5 and 7e// = 1. We find 

deff = 4 - /iH — h 

O 

«.„ = |-?^A^ (10.29) 

The model with N = 00 [49], which is related to the spherical model and the Bose 
gas, yields 

P{h, z) = -e{z)h + h'^ (10.30) 
^^2{h,z)=h (10.31) 

7^(/i,^)=0 (10.32) 

In this case the one loop expressions are exact and the effective exponents interpolate 
between those of the exact model in d and d — 1. We find 

4// = 4 - /i (10.33) 

the other exponents can then be written in terms of dgff. Specifically 

^eff = (10.34) 

7e// = (10.35) 

The exponent /?e// does not exhibit a crossover, retaining the value /?e// = \ throughout 
the crossover. 

= (10.36) 
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The specific heat effective exponent 



«e// = (10.37) 



is no longer singular for d^jj < 4 and the other effective exponents do not exhibit a crossover 
for d > 4, as the only fixed point of the model is aX h — 0. 

The results for other physically relevant systems such as the XY-model {N = 2) and 
Heisenberg model (A^ = 3) can be easily determined from the general results above. The 
latter two models are of particular interest for d — 3, but exhibit some special features 
which go beyond the scope of the present article, such as the Kosterlitz-Thouless transition 
in the former and the absence of spontaneous symmetry breaking in the latter. Our 
techniques can be simply adapted to the latter, in the guise of the non-linear a model, 
leading to a /? function that crosses over between the zero temperature fixed point as 2; — > 
and the N = 3 three dimensional critical point as 2; — > 00, further details will be presented 
elsewhere. 



We present our two loop results using Pade resummed Wilson functions in Figures 2 
through 9. In all graphs the horizontal axis is ln(^^/L), the different curves correspond to 
N = (polymers), N = 1 (Ising model), N = 2 (XY-model), N = 3 (Heisenberg model) 
and N = 00 (spherical model). The logarithmic corrections to scaling at the bulk end 
are clearly visible, their magnitude is as expected from four dimensional calculations. All 
curves are with the boundary condition /i = 1 at ln(^^/L) = —20, the value of h at the 
initial scale parameterizes different possible crossover curves but all curves asymptote to 
the same form. In Figure 4 we plot ?7e//, the exponent which governs the fall off in critical 
correlations at T = Tc{L). This exponent is not a monotonic function of N but attains 
a maximum for some value of between N = —2 and N — (yo where it is identically zero. 
This is the least accurate of our exponents and the peak appears to be at iV = 1, though 
more accurate values for this exponent suggest it occurs at higher values, probably = 3. 
Figure 7 shows 5e//, note that excepting the cases N — —2 and N — 00 the curves are very 
robust to changes of N. The crossover in the cases N = —2 and N = 00 arises purely due 
to the change in defj. Once again the curves are not montonic functions of there being 
a minimum in the vicinity of N = 3. Figure 8 shows a plot of the effective specific heat 
exponent aeff which measures how the singular part of the free energy changes as F or T 
varies. The extra case = — 2 is added here, since, in the case of dimensional crossover 
it is distinguishable from the Gaussian model due to the fact that 7;^ for the latter is zero 
whereas for the former it is non-zero being a measure of the changing effect of the leading 
irrelevant operator. Not only does one see the change in sign of the specific heat exponent 
as a function of A" but one also sees that the effective specific heat exponent can change 
sign as a function of ii/L. This is quite pronounced for the XY model which starts off 
positive, increases then turns negative at ~ 100 It would be interesting, based on the 
Harris criterion for the relevance or irrelevance of weak disorder, to see whether disorder 
could change from being irrelevant to relevant as a function of size. In Figure 9 we plot 
^\ — 4 — deff which also gives information about the effective dimensionality of the system. 
Notice that 7;^ is very robust to changes in A^, varying very little across the entire range 
of A'", [—2, 00]. The asymptotic values of the effective critical exponents and associated 
quantities in the three dimensional regime are tabulated below. 
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Asymptotic Critical Exponents 


N 






h 


7e// 


^eff 




- 2 


0* 


0* 


1.800 


1* 


0.5* 


0.5* 


- 1 


0.0200 


0.145 


1.820 


1.088 


0.550 


0.351 





0.0295 


0.277 


1.785 


1.175 


0.596 


0.211 


1 


0.0329 


0.388 


1.732 


1.257 


0.639 


0.083 


2 


0.0332 


0.479 


1.675 


1.330 


0.676 


- 0.029 


3 


0.0322 


0.552 


1.621 


1.395 


0.709 


- 0.126 


4 


0.0305 


0.611 


1.573 


1.451 


0.737 


- 0.211 


oo 


0* 


1* 


1* 


2* 


1* 


-1* 



* These values are exact. 



All these values are in very good agreement with corresponding high temperature series 
and experimental results. We believe the entire crossover curves are of similar accuracy. 



§§ 10.4 The Effect of Boundary Conditions on the Crossover. 

Let us now consider the case of other boundary conditions, more detailed results will be 
presented in a separate paper. Different conditions constitute an environmental probe 
which can be used to differentiate between systems as they fall into different crossover uni- 
versality classes. We begin with Dirichlet boundary conditions, but will only be interested 
in "bulk" not surface physics. We will confine our considerations to one loop for simplicity. 
In the three dimensional layered system we are again in the fortunate position of being 
able to present the one loop expressions in terms of elementary functions. The floating 
coupling is 

where y — {z^ — 7r^)2. It is necessary to consider two branches of the function y{z)^ 
however, even though y{z) has a branch point h{z) is analytic in z. The floating fixed 
point is simply h* = e{y) where 

e{y) = 1 + ^ + 2(1 + ^)-^ ^ (10.39) 

We will leave the reader to derive the corresponding one loop expressions for the effective 
exponents in terms of the floating fixed point. Clearly they will interpolate between exactly 
the same asymptotic values as for the fioating coupling derived from integrating the P 
function above. 

In the case of anti-periodic boundary conditions the extremal values of these functions 
are the same, they do however differ in the crossover region. Explicitly for d = 3, to one 

loop, one finds h{y) = {1 + — ). The floating flxed point is once again £(0) 

where this time 

3^_(feW)ta^ (10.40) 

y^ (smhy-y) 

Note that even though y again goes through a branch point, as z ranges between and oo, 
both (10.39) and (10.40) remain single valued. As was the case with Dirichlet boundary 
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conditions the effective exponents crossover between the same asymptotic values, however, 
the effective exponents between the asymptotic values are different. 

At the one loop level the expressions for general d, for Dirichlet boundary conditions, 

are 

eiz) = 5 - d - (7 - d)^ (10.41) 

7r2(n2-l)- 



n=l 



and for the universal one loop floating coupling 



00 (d-7) 

hiz) = (5 - rf) ^ ^ (10.42) 



n=l 

For antiperiodic boundary conditions one finds 



d-9 



e{z) = 5-d-{7- d)^^ (10.43) 



n=— oo 



^2 



and finally for the fioating coupling 



00 (d-7) 

7r^n(n+l)N 2 

;2 — J 



/^(-) = (5 - ^ (10.44) 



n=— oo 



Plots comparing the effects of different boundary conditions for a three dimensional 
layered Ising system can be found in graphs 12-16. The main interesting feature is in the 
case of Dirichlet boundary conditions, where there is a "dip" in the curves as one approaches 
the three dimensional regime. Thus the curves arc not monotonic, the effective exponents 
approaching their three dimensional asymptotic values from below. The origin of this effect 
lies in a similar characteristic dip present in the behaviour of £, and is a manifestation of 
the power law decay of the effects of the boundary in the case of Dirichlet boundary 
conditions as compared to the exponential decay for periodic boundary conditions. In the 

case o{ d = A the dip in e persists and indicates that for ^ ^ 1 the effective dimensionality 

is actually above the upper critical dimension! This has some very interesting ramifications 
which we will return to at a later date. In the case of antiperiodic boundary conditions 
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the effect is in the opposite direction to that of Dirichlet. For d = 3 the crossover to 
two dimensions is completed at larger values of L/^i when compared to either Dirichlet 
or periodic boundary conditions. At the three dimensional end there are persistent tails 
relative to periodic boundary conditions. Once again these are due to the power law as 
opposed to exponential decay of the boundary effect. For d = 4, in distinction to the case 
of Dirichlet boundary conditions, de// never exceeds the upper critical dimension. 

§§ 10.5 The Uniaxial Dipolar Ferromagnet 

Another model that falls into the same class as the considerations of this section is that of 
a uniaxial dipolar ferromagnet [19] [20]. In this case the propagator takes the form 

+ ao^ (10.45) 

Here the environmental factor that is followed in the crossover is ctg, the coupling strength 
of the long range dipole-dipole interaction between the spins. Once again the formal 

structure of the diagramatic series is the same. The diagrams differ in their dependence 

1 

on = CaoQ^o ^0^5 hence s{z), fi{z) and f2{z) are different functions, interpolating 
however, in the ai = 0, Q!2 = Pc'*^ prescription, between the values 1, | and ^ respectively 
for z = 0,to finite values close to those of the four-dimensional problem treated above. The 
numbers are different however, reflecting the absence of some angular integration variables. 
We have not as yet examined the two loop case in detail and reserve further comment till 
a future date. Let us therefore restrict our considerations to one loop. 

We find that ^ 

The solution of the running coupling equation is given by 

h = y ^ (10.47) 



There is an interesting distinction between this situation and that of the dimensional 
crossover, in that in the case of the uniaxial problem there is a universal form for the 
crossover from the three dimensional fixed point to the pseudo-four dimensional fixed point. 
The reason for this is that the true microscopic structure of the model is three dimensional, 
and therefore the logarithmic growth into the UV does not continue ad infinitum. The 
universal part of the crossover curve can be obtained as in the case of three dimensional 
problems by taking the limit — > oo. We find this universal coupling is given by 

h(z) = ^ (10.48) 

VrT^ln(i+4+i^) 

We present a comparison of j for the four dimensional dimensional crossover in a layered 
geometry with periodic boundary conditions with that of the uniaxial crossover in figure 
23. The initial coupling for the layered geometry is /ig = 1 for ln(l/2o) = ~20 and the 
curve for the uniaxial problem is the universal crossover exponent derived from (10.48). 
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In both cases we plot i^e// = ^zry-j, where 7^2 is retained to order h without further 

expansion of the denominator. We choose this form since it is closer to the two loop 
results as demonstrated in Figure 10. Figure 21 clearly shows that the uniaxial crossover 
and that of the layered geometry are in different crossover universality classes. 

As shown by Kamien and Nelson [50] recently, critical phenomena of, single compo- 
nent and binary mixture in directed polymers in an external field, strands of dipoles in 
both ferro- and electrorhealogical fluids, and flux lines in high-T^ superconductors, arc 
all in the same universality class as three dimensional uniaxial dipolar ferromagnets and 
ferroelectrics. It would be interesting to compute the effective critical exponents for the 
crossover behaviour in these systems based on the methodology here. 



§ 11. CROSSOVER IN THE QUANTAL ISING MODEL 

In this section we wish to consider yet another environmental factor — hi The real world 
is quantum mechanical, although this is not always a dominant feature governing the be- 
haviour of systems, the macroscopic world typically being governed by classical laws. How- 
ever, at low temperatures quantum effects play an important role, and typically quantum 
fluctuations need to be taken into account. When one is dealing with flnite temperature 
systems there are also thermal fluctuations. An interesting question therefore, especially in 
the theory of phase transitions, is how does one crossover between the two types of fluctu- 
ations in a realistic system. For very low temperatures one expects quantum fluctuations 
to dominate, and at higher temperatures classical thermal fluctuations. 

The explicit model we consider is an Ising ferromagnet in a transverse magnetic held 
r, described by the lattice Hamiltonian 

^=-^E^^i^f^|-rE^f (11-1) 

ij i 

where Sf are the components of spin in the a direction at lattice site i. The spins Sf and 
Sj do not commute and it is this non-commutation that gives the additional complications 

of the model over that of the usual Ising model. The model is known to have a second 
order transition, at a temperature T'c(r) for F < Fc, while for F > Fc no ordered phase 
occurs at any temperature. 

As we are only interested in the critical properties of this model here we will use the 
LGW Hamiltonian derived by Hertz [43] using the Stratonovich-Hubbard transformation, 
but now a time ordered version where time ordering is with respect to imaginary times 
between and p. The explicit Hamiltonian is 

XT " (11-2) 



where we have dropped operators irrelevant in the critical region, + = T — Fc, 
where Fc is the critical transverse field in the mean field approximation. The sum is over 

the Matsubara frequencies and T = where a = 4^Qijr tanh(/3r) ^ inverse 
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temperature. We have also used Kq — Kq — ^o^q^ + 0{q^) for the Fourier transform of 
the Ising spin-spin couphng. Observe that a is dimensionless, once /? and a are assigned 
the units of length, since from the original Hamiltonian we deduce that J, and F have 
the dimensions of inverse length. 

In this problem the anisotropy parameter, or environmental variable, is T, which is 
essentially just the thermodynamic temperature. Note that it is entering here in a very 
different fashion to the dimensional crossover considered in section 10. There the temper- 
ature entered in the standard way in the quadratic term of the LGW Hamiltonian, i.e. as 
T — Tc{L). We now wish to renormalize in an environmentally friendly manner. We use the 
conditions (5.2 — 5.5). The condition (5.2) implies, with the multiplicative renormalization 
of g, that g is proportional to F — Fc(T), the deviation from the critical line. Note that in 
the renormalization used here T does not renormalize, i.e. it is a non-linear scaling field. 
Our task now is to calculate the solution of the f3 function and to subsequently compute 
the effective exponents using the "black box" expressions of sections 5 to 9. The only 
thing we need to remember is that the environmentally friendly propagator that enters 
these expressions is G(k, uj) = 0^2_^lj'2_^g^ where co = 2'KnT. 

Consider first a three dimensional model. From (7.1) with N = 1, and using the above 
propagator in the diagrams, one finds a floating flxed point that interpolates between the 
classical and quantum fixed points. The quantum fixed point in this case is found to be 
zero, but there exist important logarithmic corrections to scaling which are captured by the 
floating coupling and are visible in Figures 16-20. An examination of deff reveals that it 
interpolates between four and three. This implies that the effects of quantum fluctuations 
are such as to increase the effective dimensionality of the system in the low temperature 
regime. The classical flxed point is h{g = 0) = 1.732. The effective exponents can be found 
from the floating coupling expressions of section 9. Using the two loop Fade resummed 
coupling and Fade resummed 7^2 for this model, the exponents i^eff^ Vef/^ leff and a^ff 
are plotted in Figures 16-19. Once again they are in excellent agreement at the three 
dimensional end with the results of [48]. We flnd for instance v{g = 0) = 0.639. The 
other scaling exponents can be found from the scaling exponent laws. Needless to say the 
expressions found agree precisely with those obtained from a direct calculation. Using the 
expressions in the appendices the explicit two loop Fade resummed expressions for a two 
dimensional Ising model in a transverse held can be found. Note that here we are deriving 
the effective exponents at flxed temperature, i.e. we are approaching the critical line from 
the g direction. We could equally well have looked at exponents defined with respect to the 
"environmental" T direction as discussed in section 5. We will return to this in a future 
paper. 

As far as scaling is concerned the relevant correlation length here for scaling purposes is 
the physical correlation length of the system. If we had used a temperature independent 
renormalization the relevant correlation length would have been the zero temperature 
quantum correlation length. One can write the scaling functions in universal scaling form 

FW=^2^¥-^-<^^(^)(Ter) 
In particular for the susceptibility = C^^-^^^H^Ct)- In terms of effective exponents 

X"' = exp (^j' v^ff{2 - Veff)Y^ (11.3) 

Strictly speaking the above functions are only truly universal for (i < 3. For d — 3 the 
logarithmic corrections to scaling evident in the effective exponents preclude a true scaling 



54 



form. However, the scaling form in terms of the non-hnear scahng fields of section 6 is 
"universal" as the corrections to scaling are captured by these fields in their dependence 
on initial conditions. For fixed initial conditions all scaling functions are universal even in 
the presence of corrections to scaling. 

Turning now to the two dimensional model, at one loop we can give some nice simple 
analytic forms. The crossover coupling in terms of 2; = T^y is 

h{z) = 1 + ^ ^ 



sinh(2;) 



One finds that 7^2(2;) = ^h{z), hence the crossover efi^ective exponents are to one loop 
given by 



12 12sinh(2) 

, . 7 z 

leffiz) = - + 



6 6sinh(2;) 

The effective dimensionality de// = £(0) is given by 

z'^coth(^) 

e{z) = 1 + ^-r— P 11-4 
smh z + z 

which once again clearly indicates a change in effective dimensionality. In this case from 
two up to three in the deep quantum regime. The susceptibility scaling function at one 
loop is simply x = z^- If one wants the results parametrized in terms of = F — Fc(T) 
then one must invert the equation 

9i^) = (11-5) 

(tanh(2/2))V3 

For the two loop results we would have to invert the two loop version of this equation. 

If one compares any of the two loop or one loop expressions for effective exponents 
here with the corresponding expressions for a layered Ising model with periodic boundary 
conditions, one will find that they have exactly the same functional form. In fact, after the 
replacements g ^ t, T ^ they are identical. Hence we would conclude that the d — 1 
dimensional Ising model in a transverse field and the d dimensional layered Ising model 
with periodic boundary conditions lie in the same crossover universality class. 

Of course it is well known that there is a mapping between these models. This was 
shown by Suzuki [8] using a Trotter product formula, the extra dimension being of finite 
size I3h with periodic boundary conditions, where /? is the inverse temperature. In the above 
we tried not to emphasize this point for two reasons. First of all, so that we could make 
the link at the level of the effective exponents in the context of which we have discussed 
and defined crossover universality classes; and secondly, so that we could bring to light 
an important point about our formalism. Suppose that at a theoretical level we knew 
nothing about the quantum fixed point, but that experiment had told us that the critical 
exponents at zero temperature were very different to those at the classical fixed point. How 
would one go about finding a theoretical description of the other fixed point, and indeed 
the entire crossover? Within our environmentally friendly formalism the answer is clear: 
one starts at some "microscopic" energy scale Em, one writes down a Hamiltonian for the 
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"microscopic degrees of freedom" , including all relevant environmental parameters, 

one cranks the handle of the environmentally friendly RG and sees what comes out. If 
there exists another fixed point which is reached via the effects of a relevant environmental 
parameter, then it will be seen in the environmentally friendly RG flow. If it does not 
exist then it will not appear. The point is we do not need to know a priori about its 
properties, or even of its existence, these can all be deduced from the RG. What is equally 
important is that one can see when one is not implementing an environmentally friendly 
renormalization because perturbation theory within the context of the unfriendly group 
will break down. This is one's pointer as to the fact that one had missed some important 
environment dependence. 

So, we have been able to deduce effective critical exponents to what we believe to 
be a very high degree of accuracy for the three dimensional quantum Ising model using 
two loop Fade resummed perturbation theory. We also gave explicit results at one loop 
level and formal results at the two loop level for the two dimensional model. The results 
we gave above are of course also valid for non-integer dimensions. One can even extend 
the results given to a four dimensional quantum Ising model. Comparing with some of 
the past literature we find, in contrast to the claims of Lawrie [22] , that the Matsubara 
frequencies do in fact play an important role in determining the scaling functions for 
dimensions 3 < d < 4, thereby allowing us to distinguish this crossover from that of a 
uniaxial dipolar ferromagnet or of the crossover between the Wilson-Fisher and Gaussian 
fixed points. Obviously for a four dimensional model there is no crossover as noted in 
[51] in the context of the "spherical" quantum Ising model. All effects then are due to 
corrections to scaling, though these can be of great importance for particle physics models 
in the context of Kaluza Klein theory [52] . In distinction to [43] aswell we derive flow 
equations that are explicitly dependent on the Matsubara frequencies. 

There are other quantum statistical systems whose crossover behaviour we can com- 
pute simply by reading off the environment dependence from the LGW Hamiltonian and 
following our black box methodology. Naturally the results presented hold identically for 
quantum ferroelectrics as well as quantum ferromagnets. We could also, following the work 
of Hertz [43], consider (in his terminology) itinerant ferromagnets, dirty itinerant ferro- 
magnets and itinerant antiferromagnets. These will fall into different crossover universality 
classes since for these models the propagators are 

G(k,a;„)=k2 + m§ + ^ 

for the itinerant ferromagnet, while for the dirty itinerant ferromagnet 

G(k,a;„) = k2 + mg + -^ 

and for the itinerant anti-ferromagnet 

G(k,(jJn) = + mg + T\aJn\ 

These will all give rise to different crossover scaling functions. We would therefore consider 
the environmental dependence to be relevant and classify the corresponding crossovers to 
be in different crossover universality classes. Clearly there is a rich vein to be explored in 
applying our methods to the quantum-classical crossover. 

§ 12. CONCLUSIONS AND SPECULATIONS 

In this paper we have developed in some detail our approach to field theoretic renormal- 
ization of crossovers. As is well known the independence of the bare vertex functions of 
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the normalization point gives an equation which aUows reparametrizations of the vertex 
functions in terms of new renormahzed parameters, via which one tries to implement per- 
turbation theory. We have emphasized, however, that one should be very wary since many 
reparametrizations may not be very useful for a perturbative analysis of the problem of 
interest. Such an analysis is best implemented in terms of parameters which at a par- 
ticular scale describe the relevant effective degrees of freedom at that scale, in preference 
to degrees of freedom which are associated with some other disparate one. We therefore 
emphasized that the optimum choice of parameters is very sensitive to the environment, 
the point being that as the environment effects the fluctuations, the effective degrees of 
freedom are consequently environment dependent. 

Our goal was to be able to quantitatively describe scale changes in systems that exhib- 
ited crossover behaviour. Scale changes in general can be fruitfully examined by analyzing 
the RG equation in conjunction with the equation of dimensional analysis. Combining the 
latter with our notion of environmentally dependent renormalized parameters gave us a 
powerful RG equation using which we could follow with our reparametrizations how the 
couplings which were natural at a given scale changed with scale, and hence how the phys- 
ical correlation functions themselves changed. One might roughly think of the philosophy 
implemented as: 1) check out the environment you are in at a certain scale; 2) choose 
the appropriate parameters for that environment at that scale; 3) change scale, but only 
by an infinitessimal amount; 4) reparameterize to variables which are appropriate for the 
environment at the new scale. By so doing one obtains an RG the fixed points of which, in 
the space of couplings Ai, are diffeomorphic to the points of scale invariance of the system. 
If an RG is used which is independent of a relevant environmental parameter then this will 
not be true. Some of the points of scale invariance will be inaccessible to an environment 
independent RG used alone, i.e. these points could not be seen as fixed points generated 
by the group with environment independent Wilson functions. In order to access them one 
would have to supplement the RG by extra non-perturbative information. 

Sometimes environmental factors can be safely neglected, but the art is to uncover 
which ones. If one is interested in a system with two crossovers induced by two parameters 
gi and g2, respectively, and one is only interested in the crossover induced by gi, then one 
could renormalize in a gi dependent and g2 independent way, and still capture the desired 
physics. Only the full ((71,(72) dependent RG will capture the full crossover structure 
however. When something has been neglected that ought not to have been, in all cases 
looked at so far, the theory has very civilly informed us, via the total breakdown in 
perturbation theory, that something was seriously amiss. We regard this as a very general 
property of crossover systems. We believe that the canonical incantations: "divergences" , 
"strong coupling" and "non-perturbative" are, more often than not, all symptomatic of 
the failure to take into account the changing nature of the effective degrees of freedom of 
a crossover system. 

Within the context of an environmentally friendly RG, dependent on a generic crossover 
parameter g, we went on to examine, in a formal way to begin with, how our program 
would be implemented. We derived expressions for the renormalization constants and the 
Wilson functions to two loops, which are valid for a rather wide class of crossover problems. 
A particular crossover problem could be simply addressed by writing down the diagramatic 
components of the expressions with propagators appropriate to the crossover in question. 

Further consideration of the general properties of environmentally friendly RGs led 
us to discover that there were generalizations of the standard scaling laws relating the 
effective exponents of crossovers. These effective exponents also appeared naturally in the 
non-linear scaling fields which were the arguments of the scaling functions characteristic 
of the crossover. We indicated how one could generate scaling forms in terms of non-linear 
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scaling fields for the entire crossover by taking a linear scaling field associated with one fixed 
point, exponentiating it, and replacing the exponent with the integral of the corresponding 
effective exponent. In the scaling limit when all relevant crossover length scales are much 
bigger than the lattice spacing the effective exponents and scaling functions are universal. 
For cases such as crossover in a four dimensional layered system, scaling functions retain 
a dependence on the initial coupling. This is due to the presence of universal logarithmic 
corrections to scaling. There is then no universal crossover scaling curve but rather a 
family of such curves which can be parameterized by the possible couplings at a fixed scale. 
Despite the fact that thermodynamic functions do not exhibit "true" scaling behaviour, we 
find that the effective exponent laws are still valid, including the logarithmic corrections to 
scaling. The effective exponents as functions define different crossover universality classes. 
As in the standard non-crossover case the existence of scaling laws implies that only two 
of the effective exponents are independent. However, in contradistinction to the standard 
case an essential ingredient in our analysis was the role played by the crossover of the 
Wilson function for the 0^ operator. The fact that the leading irrelevant operator can 
change its degree of irrelevance has made it impossible (to date at any rate) to follow a 
crossover with changing upper critical dimension using e expansion techniques. In our 
formalism this presented no difficulty. 

The crossover was naturally captured in terms of our floating h coupling which in- 
terpolated in a smooth way between the asymptotic fixed points. A natural substitute 
for £, e{g^g) appeared in (3{h,g$,g), the zero of which could be expressed in terms of a 
floating flxed point, h*, which could be ordered in powers of this s(g^g). In this case float- 
ing effective exponents appeared which also obeyed the effective exponent scaling laws. 
The floating flxed point had the advantage of being deflned with respect to an algebraic 
equation rather than a differential equation, and was therefore easier to compute. The 
difference between crossover curves deflned with respect to the floating flxed point and the 
running coupling as a solution of the (3 function equation is due to a correction to scaling 
which is slowly varying and non-singular throughout the entire crossover. These correction 
to scaling factors are calculable in our formalism. One can think of the floating effective 
exponents as giving a reasonably good approximation to the true effective exponents, the 
former missing some of the flne detail of the crossover. 

We summarized the crossover in the dimension of the (/)^ operator in terms of (ie//, 
which can be thought of as being a measure of the "effective dimensionality" of the system, 
Ixidig) ~ '^~deff being then viewed as a measure of the the deviation of the system from its 
upper critical dimension, and a natural generalization of the correction to scaling exponent 
uj. In the case of flnite size crossover in dimensions d < 4 (in the absence of transients) deff 
can be regarded as a measure of the true change in dimensionality. This interpretation of 
effective dimensionality was also seen to be quite natural as deff appeared in the scaling 
laws involving the effective exponents /3e//, ^e// and a^ff in exactly the same place as d did 
in the non-crossover case. d*ff = ^(O) one loop can also be thought of as a measuring 
the change in dimensionality. 

We discussed perturbation theory in some detail showing how and why h or h* can 
be used to order perturbation theory. Expansions in either e or 1/N without environmen- 
tally friendly renormalization are both perturbatively useless for crossovers. Even with 
environmentally friendly renormalization an e is unsuitable when the upper critical dimen- 
sion can changes during the crossover. Similarly is equally unsuitable in cases where 
the number of components of the order parameter can change. The only valid choices 
for all the crossovers considered were h and h* . We showed that in the crossovers where 
these couplings could become ~ 1 that our methodology lent itself readily to resummation 
methods, giving expressions for the [2,1] Pade resummed floating fixed point and P func- 
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tion equation. We emphasized that in our RG methodology it was the Wilson functions 
that were perturbatively computed. Scaling functions were obtained by integration of the 
corresponding characteristic equations. With expressions such as exp j {h -\- h? ...) which 
result from solving the characteristic equations, it was important not to try and expand 
the exponentials. The job of the RG is to replace perturbation theory directly at the level 
of the correlation functions with perturbation theory at the level of the Wilson functions 
the subsequent integration (exponentiation) of which ensures that the information in the 
exponentials is highly non-perturbative. Just how non-pcrturbative is seen readily in the 
explicit diagrammatic expressions for the non-linear scaling fields in section 10. Relative 
to standard perturbation theory there is clearly a huge amount of information in such 
expressions. The moral was that one should include as much information as possible in 
the Wilson functions as this information will be "exponentiated" by the RG. What is left 
out will not be exponentiated and will have to be treated in direct perturbation theory. 

We then shifted attention away from a generic crossover to some specific examples. 
We examined in detail the case of dimensional crossover. The environmental variable in 
this case was the finite size of a layered geometry. We found that our environmentally 
sensitive RG Wilson functions depended on the thickness of the layers L. The immediate 
effect of the inclusion of L was that the differential generator associated with the RG 
now interpolated between two extremal forms, that appropriate to a bulk system in the 
L/^l — > oo limit, and that associated with a reduced system in the L/^i — > limit. This 
had the desirable effect of including both extremal exponents associated with the respective 
points of scale invariance, which were now fixed points of this RG, in a natural way into the 
Wilson functions. As a consequence of this, it was easy to see that the effective exponents 
interpolated between those of the extremal fixed points and furthermore to verify that the 
effective exponent laws were valid throughout the crossover region. 

In the context of specific two loop results, we faced the issue that the non-trivial 
fixed point of the j3 function disappears at this order. We circumvented the problem 
by following Parisi [40], and using a [2, 1] Fade approximant to re-sum the series. The 
resulting (3 function equation was integrated numerically, and the solution substituted into 
the resummed perturbative series for the Wilson functions. In the four to three dimensional 
crossover the logarithmic corrections to scaling at the four dimensional end were seen to 
naturally crossover to power law scaling at the three dimensional end. Asymptotic effective 
exponents are in excellent agreement with the Callan-Symanzik/fixed dimension results of 
Baker et al. [48] and also with corresponding high temperature series and experimental 
results [41]. 

A complimentary expansion which can be included in this approach is the expansion. 
This is an appropriate one in the case of dimensional crossover since it does not rely 
on the uppercritical dimension remaining fixed. This should prove useful for examining 
dimensional crossovers as long as the layer dimension remains larger than two. We have, 
therefore, included the N = oo results, in which case the results obtained interpolate 
between the exact spherical model results for four and three dimensions, again exhibiting 
the crossover from logarithmic to power law behaviour, as found by Barber and Fisher 
[13]. For polymer type systems {N — 0) we found similar good agreement with accepted 
results. For completeness we presented the results for the = —2 model the experimental 
significance of which we are, however, unaware. It had the curious property of having only 
a crossover in the anomalous dimension of the (f)^ operator, due to the changing upper 
critical dimension of the model. Perhaps some experimental model can be found in accord 
with this. We found specifically that the exponents cue//, Peff and 5e// exhibit a crossover 
even though i^eff, left' Veff retain their mean field values. 

Boundary conditions are another natural environmental factor which affect the critical 
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fluctuations of the system. We briefly compared the differences between periodic, Dirichlet 
and antiperiodic boundary conditions on the dimensional crossover, above the transition 
temperature, and sufficiently far from any boundaries. In the three cases the formal 
structure of the RG is identical. We concluded that the effect of periodic and antiperiodic 
was similar, there is an additional shift, however, of the critical temperature, in the anti- 
periodic case and the functions /l(^) and /2(^) slightly different, but have the same 

asymptotic, large argument values. In the case of Dirichlet boundary conditions, we found, 
interestingly, that there is a "dip" in the crossover curves, hence they are not monotonic. 
The origin of this effect lay in a similar characteristic dip present in the behaviour of £, and 
was a manifestation of the power law decay of the effects of the boundary in the case of 
Dirichlet boundary conditions as compared to the exponential decay for periodic boundary 

conditions. In the case of cZ = 4 the dip in e persisted and indicated that for ^ ^ 1 the 

effective dimensionality was actually above the upper critical dimension! 

The problem of crossover in uniaxial dipolar ferromagnets was similarly briefly ex- 
amined. Here the environmental factor is the dipole dipole interaction, which modifles 
the form of the propagator. The diagramatic series is the same, however the crossover 
functions have a qualitative difference which persists to the quasi-four dimensional end. 
A true universal crossover curve exists in which all dependence on the initial coupling to 
dissapears, in contrast to the dimensional crossover problem in a four dimensional layered 
system. 

We considered also an Ising model in a transverse magnetic fleld. For the three dimen- 
sional model wc derived two loop Fade resummed expressions for the effective exponents 
describing the quantum to classical crossover. For the two dimensional model we derived 
very simple expressions for the effective exponents. The fact that the effective exponents 
for the three dimensional model were identical with those for a four dimensional layered 
Ising mdel with periodic boundary conditions showed that these two models lay in the 
same crossover universality class. We compared and contrasted our results with others 
that exist in the literature and briefly discussed how our results could easily be adapted 
to other quantum/classical crossovers. 

Because of the generality of our approach and the ubiquitousness of crossover behaviour 
this work has natural extensions to many other systems. To list but a few that readily 
come to mind: one would like to incorporate the effect of surface Hamiltonians. One 
would also like to be able to include the effect of more complicated backgrounds, such 
as the effect of the earth's gravitational field in the case of binary fluids. A background 
of instantons, or quasi-particles, or vortices for more general Hamiltonians are similarly 
natural directions in which to extend the current analysis. The effect of these we would 
expect to be most naturally encapsulated in a position dependent RG (the same would 
be true for the application of the methodology to a generic curved space necessary for 
doing early universe physics). The inclusion of dynamical effects would likewise be a 
natural thing to do. One would like to see the crossover of dynamical scaling effective 
exponents as the neighbourhood of different fixed points is reached, the RG trajectory 
could naturally be linked to time evolution, an idea which could be very useful in black 
hole formation or infiation where time gets linked to scale in a very natural way. In 
general one would expect to have a direction, position, and time dependent RG which 
responded to the anisotropic and time dependent nature of a generic environment. In the 
neighborhood of a given fixed point, however, the universality of that fixed point should 
become apparent as the associated degrees of freedom begin to fiuctuate critically. There 
is also a natural rich mathematical structure apparent in our formulation, associated with 
the geometry of the space of coupling constants M., which we plan to illucidate further in a 
future publication. Another interesting area which caught our attention recently is that of 
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microphase separation [53] , which describes the critical fluctuations in diblock copolymers 
[54] and chemically crosslinked two-component networks [55]. Here the propagator of 

interest is 1/ (A;^ + + p-) where m? is the deviation from the critical value in the Flory 

interaction parameter, and C — Rq, Rq being the radius of gyration, in the case of diblock 
copolymers, and is the elastic constant in the case of crosslinked two component networks. 
In our black box methodology all one has to do to determine the effective critical exponents 
is to substitute the above propagator into the diagramatic expressions we have given. 

Our approach has natural applications in a broad range of areas, from cosmology to tur- 
bulence, where RG ideas have proved useful. The generic features we have outlined will be 
similar. The diagramatic structure is determined from the underlying Hamiltonian/action. 
The appropriate renomalization should then be addressed, taking due care to incorporate 
the environmental factors necessary for the most efficient treatment of the processes under 
study. One might be led by effective Hamiltonians which give good phenomenological de- 
scriptions of the physics in different regimes. If one is sufRcently insightful the RG chosen 
would be the one which tracked the trajectory between the microscopic theory and the 
effective theory. Being so insightful is of course frequently impossible. The beauty of our 
approach is that one does not in principle need to know the other end. The response of 
the theory to its environment will tell you if you monitor it sufRcently closely. This in 
principle could be of help in theories where one does not know the appropriate degrees of 
freedom at one end of the crossover such as quantum gravity. 
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Appendix A. EVALUATING TWO LOOP DIAGRAMS IN A d DIMEN- 
SIONAL LAYERED GEOMETRY 

The task of this section is to evaluate the diagrams entering the RG equations in the 
dimensional crossover problem. We first list some basic formulae that prove useful. The 
introduction of Feynman parameters in combining two propagators amounts to 

1 _ r(ai + a2) f\ x-^-\l-xr^-^ 
Dl^Dl- r(ai)r(a2) k {xD^ + (1 - x)D2r+-- ^ ^ ^ 

A second useful identity is 

dJ'q 1 _ r(z/-§) 1 



(27r)«[g2 + ^2]- r(zy)(47rr/2 ^--f 



(A.2) 



We begin by evaluating the diagram O- Our diagramatic notation here is that the 
dots on a diagram represent the location of insertions in that diagram. The insertions 

could arise from different sources. O foi" instance arises at 0(A) in F^^'l) and at 0{}?) in 

r(4). 

= \y I r-, — ^ rv (A-3) 
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where 

and k can be composed of several external momenta depending on the form of the insertion. 

For periodic boundary conditions uj{n) = (27m) , n = 0, ±1, ±2 . . .. We consider only n = 
on the external lines. On introducing a Feynman parameter, performing the integration 
over q, and taking advantage of the translation invariance of the range of this integral, O 
becomes ^ ^ 

^ ^ " [a;(l - xj/c^ + T + 

For our generic normalization point with symmetric point momentum 

Pi -Pj = ^(4(^ii - (A.6) 

(the Pi are the external momenta entering the vertices which combine to form k), and mass 
point 

r = a2 (A.7) 

we have 

Evaluating the action of A^_^ = k-^ — (4 — d) we get 

For the symmetric point normalization conditions, and zero mass, we have ai = 1 and 
a2 — 0. Thus the standard symmetric point prescription yields 

r(5-d\ rl / X (d-5)/2 

and 

For a mass normalization condition at zero symmetric point momentum, we are working 
at = and a2 = 1, in which case O becomes 
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and 



The advantage of the mass point prescription is that there is one less Feynman parameter 
integral to do. 

We see that A(^4_^-jO is in fact finite when d < 5, and from the considerations of 
section 4 determines the 1-loop RG equations. A useful observation is that the asymptotic 
limits of all of these expressions for large and small kL can be obtained very simply from 
one another. Observe first, that in the limit kL 0, only the n — term contributes, 
thus this limit is obtained by retaining only the n = term in the above expressions. In 
contrast in the limit kL ^ oo the sum over n becomes an integral, which can be thought 
of as another momentum integral times a factor of kL. Thus wc can account for this limit 
by multiplying the expression in the limit kL = (obtained by retaining only the n = 
term) by k,L and replacing d — 1 by d. This is equivalent to viewing the kL — > limit of 
the d dimensional layered geometry as the limit kL ^ oo of a, d — 1 dimensional layered 
geometry. The limiting expressions don't know the difference. Of course, one can obtain 
the asymptotic expressions by direct methods if so desired, which allows the computation 
of leading asymptotic corrections. For our purposes the expressions themselves and their 
asymptotic limits are suflacient. Using the above reasoning we see that the asymptotic 
limits are 

for kL — > and 

for kL — i> oo. 

In the RG equations the quantity 



^(4-d)0(«^) = ~2 3/2 dx[x{l - x)ai + a2\ 



(d-4)/2 



s{kL) =4:-d-K-^ H^A-dO] (A. 14) 

dK 

is of interest. This expression is clearly finite and, by our rules for relating the asymptotic 
values, interpolates between A — d and 3 — d as kL ranges from to oo. 

Next we turn to the evaluation of the "cone" diagram {[), which appears in the two 
loop renormalization of F^^) and F^^'-*^), and is given by 



X 



(-LY „ti. J'n U 1,1 + r + + r + ^] ^^^^^ 

l(ki - + r + Sgi^l l(k2 + ,1 - ,2)2 + r + ^*^] 



where ni2 = ni — n2, and the diagram is evaluated with zero external discrete momenta. 
Note that the diagram only depends on the two momentum variables ki and k2 since we 
must have overall momentum conservation and there are only three vertices. 
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We proceed to evaluate this diagram by first combining the second and last propagators, 
performing the change of variables q'2 92 + x{k2 + qi), to obtain 



r(^) /■! , 1 



(A.16) 



dx- 



U^){d~l)/2 Jq ""^ a;(ni,n2,a:) i(5-^)/2 



where U!{ni,n2,x) = a;(ni2)a; + a;(ni)(l — x). Next, grouping the other two propagators 
we obtain 

1 



Wi + r + ^][{ki-qir + r + ^ 



[ dz 2 



(A. 17) 



Combining (A.16) and (A. 17), translating the qi integrand, and introducing a third Feyn- 
man parameter y to combine the two remaining propagators, yields 



r(V) 



Ad- 



[x{l-x)f-^y\l-y)y(^-d)/^ 



ni,n2 ■ 



[(1 - y)ql + y{k3 + zki + qif + z{l - z)kl + tt{x, y) + ^^^^^ 

where r(x, y) = —p. ^ 

performing the integration over q\ we have 



a;(ni,n2,a:,t/) |(^ 

(A.18) 



^ (1 - y) and Lj{ni,n2,x,y) = + a;(ni)(l - y). Finally, 



4 

ni,n2 



(A.19) 



where 



[a;(l-a;)]('^-5)/2(l-y)y(3-'^)/2 



ni,n2 



(1 - y){y{zki + k2Y + z{l - z)kl} + tt{x, y) + ^^^iM^ 



uj{ni,n2,x,y) 



(5^ (A.20) 



At our generic normalization point, kf = ai, k^ = |ci;i, kik2 = —^Oii and r = q;2; thus 
we find 



r(5 - d) 
(«L)2(47r)('^-l) 



ni,n2 



T) ^ da; / / d2; 



^0 



[a;(l - a;)]('^-5)/2(l - y)yi^-d)/2 
's{y,z)a, + r{x,y)a2 + ^^^^^^^ii^ 



(5-d) 



(A.21) 
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where s{y, z) = z{l — z){l — ip') + — y). This gives 



A 



2(4- 



r(6 -d) ^ 
(KL)2(47r)(^-l) ^ 



ni,n2 



4 

ni,n2 



(A.22) 



where 



,4 _ [x{l - - + r(x, y)a2} 



(A.23) 



For a mass independent symmetric point normalization one simply sets ai = 1 and a2 = 0, 
and for standard mass normalization at zero symmetric point momentum one sets cui = 
and 0:2 = 1 

The expressions for ^2(4-(i)'^ ^^^^^ have a simple pole at d = 4 arising from the second 
order pole in the preceding expression. The pole is cancelled, however, by the diagram 
which also appears in the two loop renormalization of T^^\ The expression for the 
combination ^2(4-d)(4) ^ ^O^) is non singular in the limit d ^ A. It is clearly dcsircablc 
to combine our expressions so as to obtain manifestly finite quantities. To this end we 
note that 5^2(4-d)0^ — 0^(4-(i)0 'the expression 



0^(4-d)0 = 



2r(6 - d) 



y2 f f dy f 
77t^ Jo Jo Jo 



(nLfiAnf-'K'^Jo Jo Jo 



dzip^ 



ni,n2 



(A.24) 



where 



[x{l - x)f-^y^y(^-dy\l - yf-^y^{z{l - z)ai + r{x,y)a2} 



ni,n2 



(6-d) 



(A.25) 



where s{x, y) = U + z{l — z){l — y) and uj{ni, 72-2, y) 



x{l—x) 



+ u{ni){l- y). 



The two expressions for ^2(4-(i)4) 0^(4-d)0 ^"^^ sufficently similar that we 
can combine them under the Feynman parameter integrals. The resulting expression is 



^2(4-d)(^-^0') 



2r(6 -d) ^ 
(KL)2(47r)^'^~^^ ni,n2 



dx I dy I dz^ 



-ni,n2 



(A.26) 



where 



ni,n2 



ni,n2 v^ni,n2 L'^'V-^ •*'/J y ^ 

(1 - + m2] (1 - - z)ai + ^2] 



sai ra2 -\ -j^zj^ 



(6-d) 



sai -\- ra2 i ^^2^2 



(6-d) 



(A.27) 
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The singularities of the integrands near y — each give rise to pole contributions in the 
individual diagrams, but these now cancel at the level of the combined integrand and the 
expression is regular at d = 4. 

Now, by our rules for extracting the asymptotic values of these expressions we obtain 



^2(4-rf)(^-^O^)UL=0 



2r(6 - d) 



(47r) 



(d-l) 



dx 



dy / dz^ofi 



(A.28) 



where 



"(1 - y)[sai + ra2] (1 - y)^^-'^^/^[z{l - z)ai + 02] 



(A.29) 



and the kL — > 00 limit is obtained by replacing d — 1 hy d and multiplying by k'^L'^ (one 
power of kL for each sum). Note that 



fl 



(A.30) 



is a finite expression for all values of kL and interpolates between ^ and ^ as kL ranges 
from to cxD for d — 4, ai — and 012 = 1, and more generally it interpolates between the 
values for cZ — 1 and d as kL ranges from to 00. 

Let us now consider the diagram arising due to wavefunction renormalization. 



^'L' rS. L L [ql + r+ ^] + r + ^] [{k + ,1 - q^f + r + ^] 

(A.31) 

and again ni2 = ni — n2- Introducing a Feynman parameter to combine the last two 
propagators and performing the integral over q2 we get 



(«L)2(47r)(^-l)/2 
1 



X 



dx[x{l — x){k -\- qi) +r 



u[ni,n2,x, 



{5-d)/2 



(A.32) 



Introducing an additional Feynman parameter and performing the integration over qi 
yields 



r(4 - d) 



{KLnA^)^'^-^) Jo 



[a;(l-a;)](^-5)/2yM/2 



(A.33) 
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We next differentiate with respect to k'^ and evaluate at our generic normalization point 
obtaining 



(«L)2(47r)(^-l) ,,,, h Jo [^(1 _ y^^^ ^ ^(^^ y^^^ ^ .in^n^^,y) f-d) 

(A.34) 

Finally we have 
^2(4-(i)0 = 

(A.35) 

The expression for 

. ffc^ (A.36) 

is therefore finite for all values of kL. In particular in the prescription ai = 0, 0:2 = 1 with 
d = 4 we have /2 ranging from ^ to ^ as kL goes from to 00. 

Appendix B. EVALUATION OF THE DIAGRAMS ON i?^ x WITH MASS 
RENORMALIZATION POINT 

In this appendix, working strictly on R"^ x wc give an alternative calculation of the 
diagrams O) {)) ^ ) ^ind the resulting expressions derived from them required for the RG 
equations. We calculate our expressions taking advantage of the fact that the propagator 
in three dimensions has the simple form 

G{x,y)= ^ . ^ e-H^-^l (B.l) 
4:7r\x — y\ 

In our four dimensional problem, we Fourier transform on the periodic coordinate (the 
of i?3 X to obtain the propagator 



G{x,y;n) = — ; -e V «^ (B.2) 

4:7r\x — y\ 

We will evaluate the diagrams at zero external momentum and also extract a scale k 
from each of our coordinates to work with dimensionless diagrams. Finally we observe 
that the corresponding three dimensional diagrams are UV finite, hence the diagrams are 
convergent for fixed values of labeling the discrete momenta propagating in the loops. 
The four dimensional UV divergences arise in this case from the sum over rij. We will 
therefore regulate the diagrams by treating the sums not as being sums over all integers 
but as sums ranging from —V to V where N is some very large but finite cut off value. 
The combinations of these diagrams that enter the RG equations are strictly finite and so 
no restriction on the sums in this case is necessary. 
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We begin with the one loop contribution 



. (B.3) 



n 



J {4nf 



where the translation invariance of the propagator is used to cancel the volume of i?^, 
which we formally took to be V. Transforming to spherical coordinates this integral can 
be evaluated. We find 

1 1 



= — r (B.4) 



Therefore, setting r = 1, i.e. a2 = 1, and acting with Aq we get 

The restriction on the sum over n can be removed for this contribution since it does not 
lead to any divergence. 

We next need to evaluate the diagram {), which is given in position space by 
i) = ^Yl J^^^J dx2G{xi,X2;m2 + mi2fG{x2,xs;mi)G{xi,xs;mi) 



ni,n2 



J (fixi J (fix2 J (fix^ 



:,-{m2+mi2)\xi-X2\ p-mi\x2-X3\ p-mi\xi-X3\ 



(B.6) 



^(^Trr^^^y J J \x1-x2V \x2-xz\ \xi-X2\ 

where we have introduced the notation 

uin) 

with ni2 = ni + 712- Again the sums are treated as being cut off at n = A?". 

This diagram, which has the geometrical shape of a cone, can be evaluated most simply 
by making a judicious choice of coordinates, taking full advantage of the translational 
invariance of the propagators. The most convenient is to choose the origin to be at one 
of the vertices, not at the apex of the cone. We choose the origin at X2 = 0. With this 
choice the integration over the location of this vertex cancels the volume factor arising in 
the definition of the vertex. The cone as a geometrical object in this setting is in fact a 
triangle since two of the propagators coincide and as such is specified by an angle and the 
length of the two sides. Thus if the angle between the points xi and ^3 is 9 the graph 
becomes 



E-i j-co t-co / , X p-'mi\xi-X3\ 
/ dxi dx3X3e-(^2+^^2)^ie-^i^3 / ^ ^ (b.7) 
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where now — xal = Wa;^ + — 2xix^iJ, The integral over /x = cos 6* can be done using 



djj, 



-m^ {a—h\i) 



bm 



-■m\/a+b\ 



_1 ^fa^^hjl 



(B.8) 



Going to the coordinates x = ^{xi + x^) and z = ^l+xl expression for the diagram 
becomes 



1 



dx I dz- 



1 „-[{m2+mi2){l+z)+mi{l-z)]x 



1 + z 



The integral over x can now be performed to get 



-2mi\z\x 



— e 



-2miX] 



(B.9) 



E 



1 



ni,n2 



1 



dz 



1 



1 



(47r)^mi 7-1 (1 + -s) yM + {M -2mi)z + 2mi\z\ M + 2mi + {M - 2mi) z 

(B.IO) 

where M = mi + m2 + ?Tti2- This integral can be performed by dividing the range of 
integration up and noting that 

dz _ 1 



JO 



I 
I 



dz 



{l-z){a + bz) a + b 
1 1 



(1 + z){c + dz) c-d 




Care is needed with the limits as there are apparent logarithmic singularities. These 
however cancel. Completing the evaluation of the diagram one obtains 



^ ~ 2(47r)2 ^ 



ni,n2 



miM 



(B.ll) 



Note that if we keep our cutoff N fixed, then in the limit kL — > the expression reduces 
reduces for r = 1 to {) = l/(967r^). Evaluating Aq^ as with AqO by setting r = 1 and 
then acting with Aq we obtain 



ni,n2 ^ 1 



(B.12) 



where 



Similarly 



1 _ 1 ^ 1 ^ 1 



M mi 1712 fni2 
1 



(B.13) 
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We are now in a position to put the two diagrams together to evaluate Aq{^ — ^O^)- We 
first evaluate 



1 -.2^ 



1 



E 



J_ 1 1 



(B.14) 



where there is no need to truncate the sum since this combination of diagrams is no longer 
UV divergent in four dimensions. 

Putting the expressions together, we obtain 



fi = 



1 



■) + 



1 



J_(J_ _^ 

mf M 2m2' miM^M 



n 



(B.15) 



The final diagram to be computed is the diagram contributing to wavefunction renor- 
malization. We evaluate it in position space but with a momentum k, in the layers R^, 
flowing through the diagram. Again we Fourier transform in the layers, i.e. the S^. The 
diagram is therefore given by 

O = y -^-^ I (fize^^^G{z;mi)G{z;m2)G{z;mi2) (B.16) 
(kLY J 



ni,n2 



where translation invariance of the propagator has been used to cancel the volume factor, 
and the sums are again regulated by a large n cutoff at N. This still has a logarithmic 
divergence due to the corresponding three dimensional divergence. The quantity of interest, 
to us at any rate, is the contribution to wave function renormalization, which is obtained 
from the derivative of this with respect to k^. It is therefore this quantity we will focus on. 
Using the explicit representation for the propagator and choosing the angle between the 
in coming momentum k and the direction z so that kz = kr/j,, where fx = cos $, we have 







1 



E 



ni,n2 ^ 



oo 



dr 



-Mr 



d_ 

di^ 



djie 



ikr/j. 



Noting that 



^^,^fcrM = 2sin(fer) 
7-1 kr 

'j: rdye-fy^c^) 



we have 

This can easily be evaluated using integration by parts yielding 

1 



(B.17) 

(B.18) 
(B.19) 



ni,n2 







327r2«2L2A;2 



E 

nin2 



M 1, , 



(B.20) 
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Noting that as A; — > 

-^ + T'^" (m) = -3(m) +■■■ 

Evaluating this in the hmit k ^ yields 

^ ^ ~ ^ 967r2«2L2M2 ^^-^^^ 

nin2 

Finally the contribution needed for the RG equations is 

ni,n2 

The interpolating function /2 — 7^-, remembering that here d = 4 is then given 

by 

4 M^M 

h = (B.23) 

( m'^ 

\ n / 

With T = 1, in the limit kL 0, only the ?ij = modes contribute, therefore m ^ 1 
M ^ 3 and M ^ |. Thus we obtain /2(0) = In the other limit, kL — > (X), the 
summations become integrals and we find /2(oo) = |. 
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FIGURE CAPTIONS 

Figure 1: Graph of £, /i and /2 for d = 4 and periodic boundary conditions. 

Figure 2: Two loop [2,1] Fade resummed floating coupling, h, ior d — 4 and periodic 
boundary conditions. N = (polymers), N = 1 (Ising model), N = 2 (XY-model), 
N = 3 (Heisenberg model), N = oo (Spherical model) and N = —2 are shown. 

Figure 3: Graph of two loop [2,1] Fade resummed u^jf for d = 4 and periodic boundary 
conditions, N = 0,1, 2, 3, oo. For N = —2 we have i'{z) = 1/2 identically. 

Figure 4: Graph of two loop [2,1] Fade resummed rjeff for (i = 4 and periodic boundary 
conditions, N = 0,1, 2, 3. For N = —2 and N = oo, r]{z) = identically. 

Figure 5: Graph of two loop [2,1] Fade resummed 7e/ / for d = 4 and periodic boundary 
conditions, = 0, 1, 2, 3, oo. For = —2, 7(2;) = 1 identically. 

Figure 6: Graph of two loop [2,1] Fade resummed /3e// for d = 4 and periodic boundary 
conditions, N = —2, 0, 1, 2, 3, 00. 

Figure 7: Graph of two loop [2,1] Fade resummed 5e// for d = 4 and periodic boundary 
conditions, N = —2, 0, 1, 2, 3, 00. 

Figure 8: Graph of two loop [2,1] Fade resummed a^ff for (i = 4 and periodic boundary 
conditions, N = —2, 0, 1, 2, 3, 00. 

Figure 9: Graph of two loop [2,1] Fade resummed 7;^(= 4 — (ie//) for (i = 4 and 
periodic boundary conditions, A'" = — 2,0,l,2,3,oo. 

Figure 10: Graph comparing different approximations for Ueff for a four dimensional 
layered Heisenberg model and periodic boundary conditions. The upper graphs cor- 
respond to two loop [2,1] Fade resummed results. In the topmost curve 7^2 was 

resummed as well as (3 and gives z^(oo) = 0.709, whereas in the lower only (3 was 
resummed and v expanded to second order in h yielding 1/(00) = 0.706. The bottom 
two graphs are the corresponding one loop results. 

Figure 11: Graph comparing two loop [2,1] Fade resummed values for the floating 
coupling and the floating fixed point. 

Figure 12: Graph of u,,jf at one loop for periodic, Dirichlet and antiperiodic boundary 
conditions, d = 3 and A'' = 1. 

Figure 13: Graph of a^ff at one loop for periodic, Dirichlet and antiperiodic boundary 
conditions, d = 3 and N = 1. 

Figure 14: Graph of at one loop for periodic, Dirichlet and antiperiodic boundary 
conditions, d = 3 and = 1. 

Figure 15: Graph of dgff at one loop for periodic, Dirichlet and antiperiodic boundary 
conditions. d = 3 and N = 1. 

Figure 16: Graph of two loop [2,1] Fade resummed 7^2 and u^ff for three dimensional 
quantal Ising model. 

Figure 17: Graph of two loop [2,1] Fade resummed a^ff and rj^ff for three dimensional 
quantal Ising model. 
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Figure 18: Graph of two loop [2,1] Fade resummed 7e// and 7;^ for three dimensional 
quantal Ising model. 

Figure 19: Graph of Veff-, 7a ^^id 7^2 at one loop for two dimensional quantum Ising 
model. 

Figure 20: Graph of two loop [2,1] Fade resummed 7e// for d = 4 and periodic 
boundary conditions showing the effect of different initial conditions. 

Figure 21: Comparison of u^ff at one loop for uniaxial dipolar ferromagnet, three 
dimensional quantum Ising model and four dimensional layered Ising model with 
periodic boundary conditions. 

Figure 22: The floating coupling h for three dimensional layered Ising model including 
crossover to mean field theory. 
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